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^H , 

^ ■ ABSTRACT 

5— ( ' 

^ , Aims. A carefully selected sample of early B-type stars in OB associations and the field within the solar neighbourhood is studied 

h^ ■ comprehensively. Present-day abundances for the astrophysically most interesting chemical elements are derived in order to investigate 

f«^ ' whether a present-day cosmic abundance standard can be established. 

Methods. High-resolution and high-S/N spectra of well-studied shaip-lined early B-type stars are analysed in non-LTE. Atmospheric 

parameters are derived from the simultaneous establishment of independent indicators, from multiple ionization equilibria and the 

hydrogen Balmer lines, and they are confirmed by reproduction of the stars' global spectral energy distributions. 

■ Results. Effective temperatures are constrained to 1-2% and surface gravities to less than 15% uncertainty, along with accurate 

-^ ' rotational, micro- and macroturbulence velocities. Good agreement of the resulting spectroscopic parallaxes with those from the new 

^T^ , reduction of the Hipparcos catalogue is obtained. Absolute values for abundances of He, C, N, O, Ne, Mg, Si and Fe are determined 

^■^ . to better than 25% uncertainty. The synthetic spectra match the observations reliably over almost the entire visual spectral range. 

r^ ' Conclusions. A present-day cosmic abundance standard is established from a sample of 29 early B-type stars, indicating abundance 

O ' fluctuations of less than 10% around the mean. Our results (i) resolve the long-standing discrepancy between a chemical homogeneous 

I , gas-phase ISM and a chemically inhomogeneous young stellar component out to several hundred parsec from the Sun, (ii) facilitate 

O ■ the amount of heavy elements locked up in the interstellar dust to be constrained precisely - the results imply that carbonaceous dust 

f~i • is largely destroyed inside the Orion H ii region, unlike the silicates, and that graphite is only a minority species in interstellar dust -, 

C/3 ' (ill) show that the mixing of CNO-buming products in the course of massive star evolution follows tightly the predicted nuclear path, 

^ I (iv) provide reliable present-day reference points for anchoring Galactic chemical evolution models to observation, and (v) imply that 

the Sun has migrated outwards from the inner Galactic disk over its lifetime from a birthplace at a distance around 5-6 kpc from the 

— . ^ ' Galactic Centre; a cancellation of the effects of Galactic chemical evolution and abundance gradients leads to the similarity of solar 

►^ ' and present-day cosmic abundances in the solar neighbourhood, with a telltaling signature of the Sun's origin left in the C/O ratio. 

j:;^ ; (abridged) 

OO I Key words. Stars: abundances — Stars: early-type — Stars: fundamental parameters — Stars: Evolution — ISM: abundances — 

r^ , Galaxy: Evolution 

in . 

en ■ 

^^ 1. Introduction with high accuracy and precision from application of various 
f^ ! r^, „ . ,,.,-,• • ,, techniques to observations of very high quality. The solar pho- 
^ , The formation and evolution of galaxies, stars, interstellar gas ^ here, its chromosphere and corona can be studied remotely 
^ . and dust, planetary systems and even ife are tightly related to ^-^ spectroscopy. Particle collection techniques from space also 
the origin and evolution of the chemical elements, and therefore ^jj^^ .„.^ .^^ measurements of the chemical composition of the 
to the cosmic cycle of matter Theories that aim at explaining ^^j^^. ^-^^ ^^^ ^^ ^^j^. ^ ^^^ tj^les. Moreover, mass spec- 
H these phenomena need to be anchored to reference values for ^. ^^ Cl-chondrites facilitates the pristine composition of 
C?. chemical abundances. The most appropriate source of reference ^^e solar nebula to be accessed (except for a few volatile el- 
abundances for the chemical composition of cosmic matter is a ^^^^^^y ^^e wealth of data gave rise to the solar abundance 
topic of vivid discussion. standard, which has been subject of active research in the past 
The chemical composition of most astronomical objects is ^^^^^^^ ^^^ -^ -^ ^^^jj ^^^-^^^ ^f continuous revision and improve- 
traditionally compared to that of the Sun because it is the closest ^^^^ ^^^^^ ^„ ^^e modeling of its conve ctive outer envelope 
star and therefore its che mical abundances can be determined ^^^ ^^e radiative transfer in 3D models (e.g.|AsElund et af^O^ 
* Based on observations collected at the Centro Astronomico His- .Lodders et al ...2009;. Caff au et al...2010i and references therein), 
pano Aleman (CAHA) at Calar Alto, operated jointly by the Max- 
Planck Instimt fur Astronomie and the Instituto de Astrofisica de Conceptually, the chemical composition of the Sun consti- 
Andalucia (CSIC), proposals H2001-2.2-01 1 and H2005-2.2-016. tutes an excellent reference for numerous astrophysical studies, 
** Based on observations collected at the European Southern Obser- of e.g. low-mass stars' interiors and atmospheres, the Galactic 
vatory, Chile, ESO 074.B-0455(A). chemical evolution of 'older' star populations, or of solar twins 
*** Based on spectral data retrieved from the ELODIE archive at Ob- (with and without planets), to mention only few among many 
servatoire de Haute-Provence (OHP). other applications. In doing this, one supposes that the Sun is a 
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typical, middle-aged low-mass star. However, the concept faces 
complications as different sets of reference values for the Sun 
are discussed in the literature. Objective criteria that allow one 
to decide which is the most appropriate on a level where the de- 
tails matter do not exist, as each of them may reproduce one or 
another aspect of observations betteiQ. Because of this we will 
refer to the generic solar values henceforth as the 'solar stan- 
dard' and specify individual sources where relevant. 

There are also examples where the Sun is considered as a 
reference because of the lack of another reliable standard, de- 
spite its chemical composition may not be representative for the 
objects of study, e.g. whenever a local and present-day refer- 
ence is required. This is because the Sun can a priori not be 
assumed to provide a standard for the chemical composition at 
present-day. The reason for this is that Galactic chemical evolu- 
tion has proceeded for the 4.56 Gyrs since the formation of the 
Sun, enriching the interstellar medium (ISM) in heavy elements. 
For instance, structure and evolution models of massive stars an- 
chor their initial composition to the solar standard despite these 
young stars are formed from the present-day ISM. Similarly, the 
solar standard is used to constrain the amount of heavy elements 
locked up in the dust-phase of the present-day ISM. Galactic 
chemical evolution models also consider the solar chemical 
abundances as representative for the 'local' Galactic neighbour- 
hood of the SurQ while it was possibly formed nowhere near 
its current Galactocentric position. Moreover, Galactic chemi- 
cal evolution models anchor the present-day composition of the 
solar neighbourhood to the solar standard in aiming to repro- 
duce present-day Galactic abundance gradients. These are only 
a few among many other applications. Establishing an appropri- 
ate abundance standard for the present-day in our local Galactic 
environment would be highly valuable either to verify whether 
the assumption that the solar standard is indeed representative, 
or to constrain various theories in the astrophysical context more 
tightly. 

Present-day chemical abundances in the solar neighbour- 
hood can be accessed relatively easily from absorption line stud- 
ies of the cold/warm ISM or from emission line spectroscopy 
of the Orion nebula. In particular the H n regions are regarded 
as privileged sites for the determination of chemical abundances 
over large distances, even in other galaxies. Unfortunately, the 
composition of the neutral and the ionized ISM gas is altered by 
depletion onto dust grains. In order to access the actual chem- 
ical composition accurately, an a-priori knowledge of the com- 
position of the dust-phase is required, which, however, is lack- 
ing. Further complications for H ii region studies are the depen- 
dence of the derived abundances on the indicators employed in 
the analysifl fluctuations of the electron temperature through- 
out the nebula, and ionization correction factors, which can be 
substantial for some chemical species. 

An ideal alternative to find a reference for the chemical com- 
position of cosmic matter are normal unevolved early B-type 
stars of ~8-18M0, which can provide simultaneously tempo- 



' The contemporarily mo st prominent e xample may be the so-called 
solar composition problem dBasu & Antiai2008.^ . the mismatch of solar 
interior models based on modem (lower) solar abundances with helio- 
seismic observations, wher eas consensus was establish ed with previous 
higher metallicity values of lGrevesse & Sauvall ([199^. 

^ The extent of the solar neighbourhood has different meanings de- 
pending on the context. For instance, in studies of F and G stars it is 
typically stretched out to several tens of pc, while here we consider it to 
extend out to ~ 400 pc, similar to studies of the ISM. 

' Recombination or c ollisionally excited lines can in dicate rather dis- 
crepant results (see e.g. lSimon-Diaz & Stasinskall201lh . 



ral (present-day) and local (birth place) information on chemical 
abundances. They can be observed not only in the solar neigh- 
bourhood, but also at larger distances in the Milky Way and in 
other galaxies. Their composition is unaffected by depletion onto 
dust grains, unlike the cool/warm ISM and the Hii regions. In 
contrast to cooler and lower mass stars they are so short lived 
(typically ~ 10'' yrs) that they have no time to travel far from their 
birth place, except for the occasional runaway star. From an anal- 
ysis perspective, spectroscopic studies of early B-type stars are 
relatively simple because their photospheres are not affected by 
strong stellar winds, unlike their hotter and more luminous sib- 
lings, or by convection and chromospheres, unlike the cool stars, 
which furthermore pose more challenges because of severe line 
blending in their crowded optical spectra. The atmospheres of 
unevolved early B-type stars are well represented by classical 
hydrostatic, plane-parallel ID-models in radiative and local ther- 
modynamic equilibrium (LTE; Nieva & P rzvbillal 120071 hence- 
forth abbreviated NP07) - though deviations from LTE (non- 
LTE effects) need to be accounted for properly in line-formation 
calculations ("Przv billa et al.ll201ll) . Moreover, they are typically 
unaffected by atomic diffusion that gives rise to peculiari ties of 
elemental abundances in many mid B- to A-type stars (e.g. lSmithl 
1996.) . In slowly-rotating stars the photosphere should also be es- 
sentially unaffected by rotational mixing with CNO-cycled mat- 
ter from the stellar core (except for some fast rotators seen pole- 
on), i.e. they should retain their pristine chemical surface com- 
position throughout their main-sequence phase. 

Despite the relatively simple physics, spectral analyses of 
main sequence B-stars have turned out to provide inconclusive 
results over the past decades, i.e. large uncertainties in basic 
stellar parameters, a tendency towards a metal -poor composi- 
tion with respect to older stars like the Sun and an overall enor- 
mous range in derived elemental abundances, challenging pre- 
dictions of stellar and Galactic evolution mo dels. For a dis - 
cussion on these topics s e e e.g. t he reviews of Herrero (2003|), 
Herrero &Lennonl (12004 . iPrzvbilla (2008) and Morel (2009) . 

In order to improve the quantitative analysis of this type 
of stars we have extensively updated the spectral modelling by 
constructing robust model atoms for non-LTE line-formation 
calculations and implemented a powerful self-consistent anal- 
ysis technique, which brings numerous spectroscopic parameter 
and abundance indicato rs into agreem ent simultaneously, start- 
ing with H, He and C fNieva^OOT). First applications of our 
method to carbon abundances in a small sample of stars in the 
solar neighbourhood (Nieva & Przvbilla 2008, NP08) provided 
an unprecedented small scatter an d an absolute average abun- 
dance similar to the solar value by lAsplund et al.l (l2005l) . solv- 
ing for the first time the above mentioned problems. Our efforts, 
also on other metals (Przybilla, Nieva & Butler 2008, PNB08), 
have provided highly-promising results so far, i.e. a drastic re- 
duction of statistical and systematic errors in stellar parameters 
and chemical abundances: uncertainties as low as ~1% in effec- 
tive temperature, ~10% in surface gravity and ~ 10-20% in ele- 
mental abundances have been achieved — compared to ~5-10%, 
~25% and a factor ~2, respectively, using standard methodO- 
The initial sample of 6 stars from OB associations and the field 
in the solar neighbourhood covering a broad parameter range 



"* Note that typical uncertainties for effective temperature given in the 
literature are <3-5%, however these do not take into account the large 
discrepancies - up to ~4000 K - that exist between different tempera- 
ture determination methods, see Sect. 5.2 of NP08 for a discussion. In 
the extremes, error estimates of up to 40-60% in surface gravity and 
a factor 3-5 in elemental abundances are found in the literature (e.g. 
iHunter et al.ll2007l : lTrundle et al.ll2007h . 
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Fig. 1. Distribution of the sample stars in the 
solar neighbourhood. Upper panel: Galactic 
plane projection, the Galactic centre lies in 
direction of the bottom of the panel. Lower 
panel: rotational plane projection, the Northern 
Galactic Pole lies in direction of the top of the 
panel. The sample stars that are analysed quan- 
titatively here are marked by dots, open boxes 
mark candidate SB2 systems and open circles 
stars that were excluded from the analysis be- 
cause of other peculiarities. Numbers/letters 
coiTespond to those of our internal numbering 
in the tables. Error bars are of spectroscopic 
distances (stars 1-20, see Table |5} and of mea- 
sured parallaxes (objects A-G). The position of 
the Sun is marked by O. 



turned out to be chemically homogeneous on the 10% level, cor- 
roborating earlier findings from analyses of the ISM gas-phase 
(PNB08). Yet, the sample size was arguably small. 

In this work we therefore present a comprehensive study of a 
sample of initially 27 early B-type stars from the solar neighbor- 
hood, aiming to re-investigate our previous findings (from a sub- 
sample of only 6 stars) on a solid statistical basis. The stars span 
about 20 000 K in effective temperature and ~0.7 dex in surface 
gravities, from the zero-age main sequence (ZAMS) to close to 
the end of the core hydrogen burning phase. Seven stars from 
the original sample had to be excluded from the analysis for var- 
ious reasons. In order to compensate for this loss, nine additional 
stars from the Ori OBI asso ciation analysed in an identical way 
dNieva & Sun6n-Diazll201 ll NSll) were included in the deriva- 
tion of the proposed cosmic abundance standard. 

The paper is organised as follows. Particular attention is paid 
to the selection criteria for the star sample and of the special care 
taken to obtain observations of high quality in Sect. |2] The com- 
putations of extensive grids of synthetic spectra and the imple- 
mentation of a semi-automatic analysis technique are discussed 
in Sect. |3] which facilitates the same quality of the analyses to 
be achieved as in our previous work (Sect. |4}. A comprehensive 
characterisation of the atmospheric properties of the sample stars 
and related quantities is provided in Sect. |5] Then, the resulting 
present-day cosmic abundance standard for the solar neighbour- 
hood is discussed in Sect. |6] Our results are put in the broader 
astrophysical context in Sect. |7] concentrating on ISM science, 
and stellar and Galactic chemical evolution. Finally, a summary 
is given in Sect. [8] 



2. The star sample 

2.1. Primary target selection 

An investigation of the present-day chemical composition of 
the solar neighbourhood requires a careful selection of the 
star sample to minimise observationally-induced systematic 
bias. We took advantage of previous studies of early B-type 
stars (lGri2sbv et al.lll992l 'Gies & Lambertlll992l, lKilia nlll992l: 
Cunha & Lambert 1994; Andrievskv et al. 1999; Adelm an et al.l 
2002; Lvubim kov et al.ll2004t IMorel et al.ll2006l) to compile our 
target list. Criteria were the stars to be 

i) bright: early B-type stars near the main sequence of magnitude 
y < 6 mag are located at distances < 500 pc, with such a bright- 
ness limit facilitating high-quality spectra to be obtained with 
relative ease, 

ii) sharp lined: low (projected) rotational velocities v sin / < 
40 km s"' allow spectral line analyses to be done at highest preci- 
sion, maximising the chances to identify line blends and to place 
the continuum unambiguously, 

m) single: this is to prevent systematic errors in the analysis 
as second light from a companion distorts the ratio of line- to 
continuum-fluxes; objects in SBl systems with much fainter 
companions or individual components in a visual binary can 
therefore still qualify as targets, and 

iv) chemically inc onspicuous: th e CP phenomenon is rare among 
early B-stars (e.g. ISmithl[l996h but spectroscopically classified 
He-weak or He-strong stars need to be excluded as segregation 
processes in their atmospheres have rendered them useless to 
trace the pristine chemical composition. 

In consequence, many of the resulting targets are among the 
best-studied early B-type stars, with multiple analyses reported 
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Table 1. The complete star sample: id, spectral type, variability, OB association membership, photometry" , and observational 
details. 



# 


HD 


HR 


Name 


Sp.T" 


Variab. 


OB Assoc. 




V 
mag 


B-V 

mag 


Date 
DDMMYYYY 


A' 


T 

^ exp 
S 


S/Ns 


Feros R 


= 48 000 


























1 


36591 


1861 




BIV 




OriOBlIb 




5.339 


-0.194 


25-02-2005 


1 


120 


530 


2 


61068 


2928 


PTPup 


B2III 


/?Gep 


Field 


+ 


5.711 
0.015 


-0.176 
0.010 


25-02-2005 


1 


180 


420 


3 


63922 


3055 




BO III 




Field 


+ 


4.106 
0.005 


-0.185 
0.007 


25-02-2005 


1 


60 


470 


4 


74575 


3468 


Q-Pyx 


B 1.5 III 


/?Gep 


Field 


+ 


3.679 
0.006 


-0.183 
0.003 


27-02-2005 


2 


60 


300 


5 


122980 


5285 


;rCen 


B2V 


/?Gep 


Up. Gen Lup 


+ 


4.353 
0.007 


-0.195 
0.006 


26-02-2005 


1 


60 


480 


6 


149438 


6165 


tSco 


B0.2V 




Up. Sco 


+ 


2.825 
0.009 


-0.252 
0.007 


24-02-2005 
25-02-2005 


1 
3 


80 


810 


FocES R 


= 40000 


























7 


886 


39 


7 Peg 


B2IV 


/JGep 


Field 


+ 


2.834 
0.015 


-0.226 
0.012 


24-09-2005 


2 


180 


530 


A 


22951 


1123 


oPer 


B0.5V 




PerOB2 


+ 


4.968 
0.008 


-0.017 
0.012 


24-09-2005 


1 


500 


280 


8 


29248 


1463 


vEri 


B2III 


/?Gep 


Field 


+ 


3.930 
0.023 


-0.210 
0.009 


14-10-2005 


1 


360 


390 


9 


35299 


1781 




B1.5V 




OriOBlIa 


+ 


5.694 
0.010 


-0.210 
0.007 


15-10-2005 


2 


2400 


330 


B 


35468 


1790 


Bellatrix 


B2III 




OriOBl 


+ 


1.636 
0.007 


-0.224 
0.008 


24-09-2005 


6 


180 


520 


10 


35708 


1810 


oTau 


B2.5IV 




Gas-Tau 


+ 


4.875 
0.012 


-0.145 
0.006 


14-10-2005 


1 


600 


280 


11 


36512 


1855 


L'Ori 


BOV 


/JGep 


OriOBlIc 


+ 


4.618 
0.013 


-0.264 
0.007 


27-09-2005 


2 


720 


500 


12 


36822 


1876 


<^'Ori 


BO III 




OriOBl 


+ 


4.408 
0.006 


-0.162 
0.013 


28-09-2005 


2 


480 


420 


13 


36960 


1887 




B0.5V 




OriOBlIc 


+ 


4.785 
0.007 


-0.250 
0.010 


28-09-2005 


1 


240 


260 


C 


37023 


1896 


(9'OriD 


B0.5 Vp 




On OB Id 




6.700 


0.080 


26-09-2001 


5 


3080 


400 


14 


205021 


8238 


/JCep 


BlIV 


ySGep 


Field 


+ 


3.233 
0.014 


-0.222 
0.006 


26-09-2005 


1 


100 


310 


15 


209008 


8385 


18 Peg 


B3III 




Field 


+ 


5.995 
0.008 


-0.120 
0.014 


08-10-2005 


2 


1800 


410 


16 


216916 


8725 


EN Lac 


B2IV 


/JGep 


Lac OBI 


+ 


5.587 
0.015 


-0.144 
0.008 


26-09-2005 


1 


600 


270 


Elodie R 


= 42000 


























17 


3360 


153 


fCas 


B2IV 


SPB 


Gas-Tau 


+ 


3.661 
0.017 


-0.196 
0.006 


14-08-2003 


4 


1260 


310 


D 


11415 


542 


eCas 


B3III 


Be? 


Gas-Tau 


+ 


3.370 
0.009 


-0.155 
0.007 


12-01-2003 


2 


300 


310 


18 


16582 


779 


(5 Get 


B2IV 


/JGep 


Gas-Tau 


+ 


4.067 
0.007 


-0.219 
0.008 


13-01-2003 


5 


1250 


310 


19 


34816 


1756 


^Lep 


B0.5 IV 




OriOBl 


+ 


4.286 
0.005 


-0.273 
0.015 


25-12-1996 


1 


1200 


250 


20 


160762 


6588 


(Her 


B3IV 


SPB 


Field 


+ 


3.800 
0.000 


-0.179 
0.003 


28-05-2003 


2 


600 


390 


E 


163472 


6684 


V2052 Oph 


B2IV-V 


/JGep 


Field 


+ 


5.823 
0.015 


0.089 
0.003 


19-08-2003 


3 


6000 


290 


F 


214993 


8640 


12 Lac 


B2III 


/JGep 


Lac OBI 


+ 


5.253 
0.018 


-0.137 
0.008 


10-01-2004 


1 


2100 


250 


G 


218376 


8797 


IGas 


B0.5 IV 




Field 


+ 


4.850 
0.009 


-0.028 
0.011 


12-11-2004 


1 


2100 


290 


lMermilliodl(ll991h:lMorel 


& Maenenad (119781) for HD 36960 and HD 37023. ("'iHoffleit & 


Jaschekl (119821). 









in the literature. The sample could therefore be viewed as 'ideal' disk, with few outliers. They are highly concentrated towards the 
for the proposed purpose. disk plane, with few objects located beyond a vertical distance 
The distribution of the sample stars in the solar neighbour- z > 100 pc. Members of the two most prominent star-forming re- 
hood is visualised in Fig.[T](anticipating distance determinations gions in the solai- vicinity are included in the sample, object C is 
in Sect. |5]i. The stars deUneate Gould's Belt both in the Galactic one of the Orion Trapezium stars, objects 5 and 6 are members of 
plane projection as well as in the vertical cut through the Galactic the Sco-Cen association. The asymmetry of the star distribution 
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Table 2. lUE spectra used in this study. 



# 


HD 


SW 


Date 


LW 


Date 


3 


63922 


P09511 


13-07-1980 


R08237 


13-07-1980 


5 


122980 


P46857 


30-01-1993 


P24819 


30-01-1993 


6 


149438 


P33008 


01-03-1988 


PI 2766 


01-03-1988 


7 


886 


P43467 


25-12-1991 


P22073 


25-12-1991 


8 


29248 


P37958 


06-01-1990 






9 


35299 


PI 8005" 


18-09-1982 






11 


36512 


P08164 


04-03-1980 


R07097 


04-03-1980 


12 


36822 


P08595" 


29-03-1980 


R07338" 


29-03-1980 


13 


36960 


P30541 


16-03-1987 


P10338 


16-03-1987 


14 


205021 


P40477'' 


28-12-1990 


P19491« 


28-12-1990 


15 


209008 


P20593'' 


03-08-1983 


R16508'' 


03-08-1983 


17 


3360 


P26535 


03-08-1985 


P09140 


21-09-1986 


18 


16582 


P29814 


05-12-1986 


P09634 


05-12-1986 


19 


34816 


P08166 


04-03-1980 


R17279 


08-03-1984 


20 


160762 


P42454 


13-09-1991 


P21228 


13-09-1991 


Notes. 


"■"' High-resolution 


spectrum. 







is because our main observing program utilised telescopes on the 
northern hemisphere (the celestial equator is roughly spanned by 
the connecting line between objects E and 9). 

2.2. Observations and data reduction 

High-resolution echelle spectra at very high signal-to-noise 
(S /N) ratio - ranging from 250 up to over 800 in B - and wide 
wavelength coverage for 27 objects were obtained, either by own 
observations, or from archives. A summary of the star sample 
and observational details are given in Table[T] Object identifica- 
tions include our own numbering scheme (to facilitate easy iden- 
tification in the figures and other tables), HD and HR numbers 
and names of the stars. Spectral types are given and an indication 
of variability, basically divided in ySCephei types and slowly- 
pulsating B-stars (SPBJ3- Membership to one of the nearby OB 
associations or to the field population is indicated. Photometric 
information covers observed Johnson V-magnitudes and colors 
B - V. Observational details like the date of observation, the 
number of individual exposures per object A^, the total expo- 
sure time and the resulting S /A^-ratio in the B band of the coad- 
ded spectrum are also indicated. Our philosophy for coaddition 
of data aims at obtaining a 'snapshot' of the stellar spectrum, 
even for variable stars. Consequently, only exposures taken dur- 
ing a small fraction of the same night were coadded, with the 
one exception of a non-variable star. The analysis will thus yield 
meaningful time-invariant quantities like elemental abundances 
or fundamental stellar parameters even for variable stars, while 
atmospheric parameters like effective temperature and surface 
gravity will be valid only for the moment of observation. The 
observational material divides into three subsets. 

The first subset includes stars 1-6, for which spectra 
were obtained with Feros ( Fiberf ed Extended Range Optical 
Spectrograph, iKaufer et all Il999l) on the ESO 2.2m tele- 
scope in La Silla/Chile. Feros provides a resolving power of 
R = A/AA X 48 000, with 2.2 pixels covering a AA resolution el- 
ement. Of the entire Feros wavelength range, only the part be- 
tween ~3800 and 8000 A meets our quality criteria for detailed 
an alysis. The stars have alr eady been subject to investigation 
in iNieva & Przvbillal dlool NP06), NP07, NP08 and PNB08. 
Nevertheless, they are included here as many details of the anal- 



ysis were unreported previously. We refer the reader to NP07 for 
details of the observations and the data reduction. 

Thirteen stars (7-16, A-C) compose the second subset, 
which were obser ved with Foces (Fibr e Optics Cassegrain Ech- 
elle Spectrograph. IPfeiffer et al.l ll998) on the 2.2m telescope at 
Calar Alto/Spain. With the instrument configuration chosen a 
2-pixel resolution of R^4Q 000 was obtained, providing wave- 
length coverage from ~3900 to 9500 A. Data reduction was per- 
formed using the Foces data reduction software (Pfeifferet al] 
1998), comprising the usual steps of bad pixel and cosmic cor- 
rection, bias and dark current subtraction, removal of scattered 
light, optimum order extraction, flatfielding, wavelength calibra- 
tion using Th-Ar exposures, rectification, and merging of the 
echelle orders. A major advantage related to the Foces design is 
that the order tilt for normalisation is much more homogeneous 
than in similar spectrographs. This facilitates a more robust con- 
tinuum rectification than feasible usually, even in the case of 
broad features like the hydrogen Balmer lines, which can span 
more than one echelle order (see Korn 2002) . Finally, multiple 
exposures of individual stars were coadded and the spectra were 
shifted to rest wavelength by cross-correlation with appropriate 
synthetic spectra. 

Spectra for the third subset of stars (ei ght stars: 17-20, D-G ) 
were extracted from the Elodie archive (Moult aka et al.ll2004l) . 
Elodie was an echelle spectrograph mounted on the 1.93 m 
telescope of the Observatoire de Haute-Provence/France, cov- 
ering a wavelength range from -4000 to 6800 A at /?«42000 
(see Baranne et al. 1996). The reduced spectra from the archive 
needed to be normalised as a last step of post-processing. We 
used a spline function to carefully selected continuum points 
for this purpose, in analogy to our procedure for the Feros data. 
Unfortunately, a reliable normalisation of broad lines in the blue 
like Hy and H6 turned out difficult, such that they were not con- 
sidered for the analysis. After normalisation, the individual spec- 
tra of one object were coadded in order to obtain a combined 
spectrum of improved S IN. 

In addition to the spectra, which constitute the principal 
data for the analysis, various (spectro-)photometric data were 
adopted from the literature for constructing spectral energy 
distributions (S EDs). M ain sources of Johnson UB VRIJHK- 
phot ometry werelMemiilliod (19911) . [Mo rel & Magn enat ( 1978|) 
and iDucatil (l2002h . and some 7i//r-magnitudes (only those 
with a quality marker 'A') were adopted from th e Two Micron 
All Sky Survey (2MASS. [Slcrutskie et al.]l2006l) . Additionally, 
flux-calibrated spectra for 15 objects as observed with the 
International Ultraviolet Explorer (lUE) were extracted from the 
MAST archiveQ Only spectra acquired using a large aperture 
were considered to avoid light loss. Preference was given to 
low-dispersion data, except for a few cases where only high- 
dispersion spectra were available. These data cover the range 
from 1 150 to 1980 A for the short (SW) and from 1850 to 3290 A 
for the long wavelength (LW) range camera. Typically, both 
wavelength ranges were observed the same day. A summary of 
the individual spectra used in the present work (data ID number 
and observation date) is given in Table |2] 

The observational material comprises one of the highest 
quality for the quantitative study of elemental abundances from 
early B-type stars in the solar neighbourhood so far Obser- 
vational bias should thus be reduced to a minimum. 



^ Pulsation periods of /3 Cephei stars vary between 
while those in SPB stars vary between ~0.5 to 3 days. 



-3 to 8 hours, 



ihttp : //archive . stsci . edu/| 
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Fig. 2. Examples of prominent spectral lines 
of the newly identified candidate SB2 systems 
in our sample, HD 22951 (B0.5 V), HD 35468 
(B2III) and HD 37023 (B0.5V). For compar- 
ison the corresponding spectral lines in the 
single star HD 36822 (BO III) are shown as 
well. The emission peak in He i /i6678 of the 
Trapezium star HD 37023 in Orion is of nebu- 
lar origin. 



2.3. Stars excluded from the analysis 

Once the spectra were reduced we undertook visual inspections 
of the data in the course of the quality assessment. Four objects 
were immediately found conspicuous, while three more were 
found to be problematic in the analysis. About a quarter of the 
carefully pre-selected sample turned out as candidate binaries 
or as chemically peculiar In the following we want to briefly 
address the reasons why these cases (stars A-G in Table [T} 
needed to be removed from our sample. 

A: HD22951 The object is a so far unrecognised candidate 
SB2 system. All the stronger sharp lines in the spectrum show 
weaker and broader absorption dips as well, shifted blueward 
from the primary's features. Examples are given in Fig. where 
the corresponding lines in one of the single stars of our sample 
are also shown for comparison. The secondary of the system 
is a cooler main-sequence star, of spectral type about B2. All 
the signatures can be viewed as undetectable at the resolution 
of classification spectra. The secondary's light renders the 
spectrum unsuited for our analysis technique and the star is 
therefore excluded from further analysis. 
B: HD35468 (Bellatrix) The star was considere d a referenc e 
for the definition of spectral type B2III by Wal bornI (Il97lh . 
However, our spectrum shows that y Ori is a candidate SB2 
system composed of two very similar components, see Fig. |2] 
Th is resolves the appare nt overluminosity of the star found 
by ISchroder et al.l (|2004|) from an evaluation of its Hipparcos 
parallax. 

C: HD37023 The case resembles that of HD 22951, see Fig.|2] 
The companion of HD 37023 is of slightly earlier type than that 
of HD 22951, as its contribution to the spectrum is stronger 
Note that while HD 37023 is a known spectroscopic binary 
(see IVitrichenkol 12002. for a discussion) the direct detection 
of a companion in the spectrurrQ offers the possibility to put 
much tighter constraints on the system, such that follow-up 
observations are recommended, like for the two previous cases. 

D: HD1 1 41 5 A preliminary analysis of the star yielded several 
inconsistencies, which either relate to unrecognised problems 
with the Elodie spectrum (in particular to the continuum nor- 
malisation) or to some peculiarity of the star. A closer inspection 
of the liter ature revealed that it is listed in the catalogue of 
Be-stars of iJaschek & Egrej (Il982h . All Balmer lines and in 
particular Ha are in absorption in the available spectrum, but we 
cannot rule out low-Ievel emission filling in the inn er line wings 

i imitating normalisation problems). Moreover, iLeone et all 
1995h describe HD 11415 as CP star of the He-weak type. 



^ The presence of light from a co mpanion in the spe ctrum of 
HD 37023 was indicated recently also by ISimon-Diad ( 120101 , his foot- 
note 3), however without giving further details. 



These factors prompted us to exclude HD 11415 from the 
detailed analysis. 

E: HD 163472 A preUminary analysis of the star indicated 
inconsistencies of the solution, which were unexpected in 
vie w of the 'normal' appearance implied by the analysis of 
Mo rel et al.l (|2006). However, a detailed study by iNeiner et al.l 
(i2003h found HD 163472 to be chemically peculiar, in particular 
it is a He-strong star. Diffusion in the magnetic atmosphere 
renders this otherwise highly interesting star useless for our 
purpose. 

F: HD214993 Th e star is one of the most intensely studied 
jSCephei pulsators (iDesmet et al.ll2009 i, and references therein). 
Our preliminary modelling encounters difficulties, which may 
be resolved assuming helium-enrichment of the atmosphere. 
This requires further investigation and in consequence we 
exclude the star from the present work program. 
G: HD218376 The star shows conspicuously broad lines. 
Whether this may be interpreted as unresolved binarity - re- 
quiring two about similar components - or whether a different 
explanation needs to be found has to be decided by further in- 
vestigations. 

3. Spectrum synthesis in non-LTE 

The non-LTE line-formation computations follow the methodol- 
ogy discussed in detail in our previous studies for H and He in 
NP07, for C in NP06 and NP08, and for N, Ne, Mg, Si and Fe 
in PNB08. In brief, a non-LTE approach is employed to solve 
the restricted non-LTE problem on the basis of prescribed LTE 
atmospheres. This technique provides an efficient way to com- 
pute realistic synthetic spectra in all cases where the atmospheric 
structure is close to LTE, like for the early B-type main sequence 
stars analysed here (see NP07). The computational efforts can 
thus be focussed on robust non-LTE line-formation calculations. 
The validity of the approach was recently verified by direct com- 
parison with hydrodynamic line-blanketed non-LTE model at- 
mospheres (NSll). This approach has been equally successful 
in improving model fits beyond th e field of massive B-type stars, 
e.g. in low-mass subdwarf B-stars ([Przvbilla et al. 2006b) and B- 
type extreme helium stars (iPrzvbilla et al.l l2005 . 2006cl). 



3.1. Models and programs 

The model atm ospheres were computed with the Atlas9 code 
(iKuruczll 1993b ) which assumes plane-parallel geometry, chem- 
ical homogeneity, as well as hydrostatic, radiative and local 
thermodynamic equilibrium (LTE). Line blanketing was realised 
here by means of opacity dist ribution functions (ODFs, Kurucj 
Il993ah . Solar abundances of iGrevesse & Sauvall (1 19981) were 
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Table 3. Model atoms for non-LTE calculations. 



Ion 



Model atom 



H 

He i/ii 
C ii-iv 
Nil 
Oi/ii 
Ne i/ii 
Mgii 
Si m/iv 
Fen/m 



Przvbilla & Butler (2004) 

Przybilla (2005) 

NP06, NP08 

Przvbilla & Butler (2001) . updated" 

Przvbilla et al. (2000), Becker & Butler (1988) . updated" 

Morel & Butler (2008), updated" 

Przvbilla etal. (2001) 

Becker & Butler (1990) 

Becker (1998), Morel et al. (2006). corrected'' 



Notes. " See Table 7 for details. * See Sect.lSTT 



adopted in all computations. The model atmospheres were held 
fixed in the non-LTE calculations. Non-LTE level populations 
and model spectra were obtained with rece nt versions of Detail 
and Surface ( Giddings 1981; Butler & Giddings 1985, both up- 
dated by K. Butler). The coupled radiative transfer and statisti- 
cal equilibrium equations were solved with Detail, employing 
an ac celerated lambda iteration scheme of Rvbicki & Hummer 
(119911) . This allowed even complex ions to be treated in a re- 
alistic way. Synthetic spectra were calculated with Surface, 
using refined line-broadening theories. Continuous opacities 
due to hydrogen and helium were considered in non-LTE and 
line blocking was accounted for in LTE via Kurucz' ODFs. 
Microturbulence was considered in a consistent way through- 
out all computation steps: in the selection of appropriate ODFs 
for realising line blanketing and line blocking in the atmospheric 
structure and non-LTE level populations determination, and for 
the formal solution. 

Non-LTE level populations and the synthetic spectra of all el- 
ements were computed using our most recent model atoms listed 
in Table [3] Updates of some of the published models were car- 
ried out introducing improved oscillator strengths and collisional 
data. These models were previously tested in NP08 and PNB08 
for early B-type stars covering a wide parameter range. A prob- 
lem with the line-formation calculations for Fe in was identified 
in the course of the present work. Several high-lying energy lev- 
els were previously erroneously treated in LTE in the formal so- 
lution with Surface, despite correct non-LTE level populations 
were provided by Detail. Higher iron abundances for the hotter 
objects (e.g. by 0. 16 dex for t Sco and <0.05 dex for the majority 
of stars) result after implementing the corrected Fe model atom, 
removing a slight artificial trend with temperature found previ- 
ously. 

3.2. Grids of synthetic spectra 

For the present work a set of grids of synthetic spectra was com- 
puted with Atlas9, Detail and Surface following the same pro- 
cedure as in our previous papers. Large independent grids of 
H/He, C, N, O, Mg and Si encompass effective temperatures 
from 15 000 to 35 000 K in steps of lOOOK, surface gravities 
\ogg from 3.0 to 4.5 (cgs units) in steps of 0. 1 dex, microturbu- 
lences from to 12kms"' in steps of 2 km s"' and metal abun- 
dances within 1 dex centered on the B-star abundance values of 
PNB08 in steps of 0. 1 dex. Hydrogen and helium abundances 
are set to the values derived by PNB08. The lower limit of the 
surface gravity for each value of temperature/microturbulence 
is constrained by the convergence of Atlas9. All grids have 
been successfully tested by reproducing results from PNB08. 
For Ne and Fe (the computationally most demanding species) 
microgrids - varying abundance only - were computed per star 




Model Atoms 



Non-LTE Spectrum 
Synthesis 



Quantitative 
Analysis 







Fig. 3. Schematic diagram of our analysis method. See Sect. IJT 
for details. 



once all stellar parameters were determined with the larger pre- 
computed grids. 



4. Spectral analysis 

Our analysis method is based on the simultaneous reproduction 
of all spectroscopic indicators (Sect. 14.2b via an iterative line- 
fitting procedure aiming to derive atmospheric parameters and 
chemical abundances self-consistently. In contrast to common 
strategies in stellar spectroscopy, this analysis technique takes 
full advantage of the information encoded in the line profiles 
at different wavelength ranges simultaneously. Integrated quan- 
tities like equivalent widths W^ are not measured in this ap- 
proach. The stellar parameters primarily derived here are the ef- 
fective temperature T^ff, surface gravity logg, microturbulence 
^, (radial-tangential) macroturbulence (, projected rotational ve- 
locity V sin / and elemental abundances s(X) - log(X/H) -1-12. 



4.1. Semi-automation of fhe analysis procedure 

The basic analysis methodology was introduced and applied to 
a strategically chosen star sample covering a broad parameter 
range in Nieva (2007), NP07 and NP08. A major complication in 
that work was that both the model atom for carbon and the stel- 
lar parameters had to be constrained simultaneously. Thus, one 
global problem needed to be solved - to select among the avail- 
able atomic data the one combination that facilitated reproduc- 
tion of the observed lines throughout the entire sample equally 
well - and many individual problems - the determination of at- 
mospheric models that described the sample stars best. The only 
way to realise this was via the computation of numerous small 
grids per star for different sets of stellar parameters and atomic 
data, and interactive work to decide how to improve based on 
the quality of match between the models and observation. Over 
100 interactive iterations were needed to find a self-consistent 
solution, i.e. selection of the best atomic data and determina- 
tion of stellar parameters and carbon abundances for six stars 
covering a broad temperature range. The benefit from this time- 
consuming approach was twofold: first, a highly robust model 
atom was provided for further applications, and second, copious 
experience was gained in the identification of sources of sys- 
tematic error that helps to minimise such uncertainties in work 
thereafter This lead, e.g. to updates of model atoms for other ele- 
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Table 4. Spectroscopic indicators for for Teg and log g determi- 
nation. 



# HD Teff H He iHe ii C ii C iii C iv O i O ii Ne iNe ii Si iiiSi iv Fe iiFe iii 
lO^K 












11 36512 33.4 . 
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15 209008 15.8 . 
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Notes. Tiie boxes denote ionization equilibria. 



ments, see PNB08 and Table[3] All in all, a solid basis for highly 
accurate (reduction of systematic uncertainty to the greatest pos- 
sible extent) and precise (low statistical uncertainties) analyses 
was thus laid. 

Once the global problem of the model atom construction and 
testing is solved, one is faced with the much easier task that only 
stellar parameters and elemental abundances need to be con- 
strained. The basic analysis methodology is outlined in Fig. |3] 
The individual steps could be realised in form of small grid com- 
putations like outlined before, but this is too inefficient for inves- 
tigating larger star samples. We have therefore replaced the cal- 
culation of dedicated small grids of synthetic spectra (compris- 
ing the boxes 'Model Atmosphere', 'Model Atom' and 'Non- 
LTE Spectrum Synthesis' in Fig. |3]l by a pre-computed compre- 
hensive grid of models (in total of the order -100000 synthetic 
spectra), as described in Sect. 13. 21 The second new ingredient for 
the present work was the adaptation of a powerful fitting routine 
for the semi-automatic comparison of observed and theoretical 
spectra (the box 'Quantitative Analysis' in Fig. [3]l. Spa^ pro- 
vides the means to interpolate between model grid points for up 
to three parameters simultaneously and allows to apply instru- 
mental, rotational and (radial-tangential) macrobroadening func- 
tions to the resulting theoretical profiles. The program uses the 
downhill simplex algorithm (Nelder & Mead 1965) to minimise 
;\f^ in order to find a good fit to the observed spectrum. 

Interactive work in some decisive steps on the analysis with 
Spas paid off as much more accurate results could be obtained. 
Crucial was the selection of the appropriate spectroscopic in- 
dicators (Sect. I4.2l i for the parameter determination which may 
vary from star to star upon availability of specific spectral lines 
(depending on stellar temperature, spectrum quality and the ob- 
served wavelength coverage). All spectral lines unsuited for 
analysis because of e.g. blends, low S/N, uncorrectable normal- 
isation problems, incomplete correction of cosmics, or known 



shortcomings in the modelling needed to be excluded. Also a 
verification and, possibly, correction of the automatic contin- 
uum placement lead to a gain in precision. Every element was 
analysed independently (passing through the loop procedure in 
Fig. [3]) and some interactive iterations for fine-tuning the param- 
eter determination were needed in order to find a unique solu- 
tion that reproduces all indicators simultaneously. This facili- 
tated also to derive realistic uncertainties for the stellar parame- 
ters, as the formal errors determined by Spas (via bootstrapping) 
were unrealistically low. Instead, the standard deviations around 
the average parameter values were adopted, as derived from the 
various independent spectral indicators. Likewise, uncertainties 
of elemental abundances were determined from the line-to-line 
scatter found from the analysis of the individual features. 

Finally, it was thus possible to derive a simultaneous, self- 
consistent solution for atmospheric parameters and chemical 
abundances, and also to quantify their statistical uncertainties. 
The novel approach provides results meeting the same quality 
standard as our previous work (for test purposes and consistency 
checks we therefore included the previously analysed objects in 
the present sample, stars 1-6 in the tables). Its advantages are a 
higher degree of objectivity than 'by eye' fits and a far higher 
efficiency, hence allowing larger star samples to be analysed. 

4.2. Stellar parameter and abundance determination 

Special emphasis was given to use multiple indicators in order 
to minimise the chance of the stellar atmospheric parameters 
and chemical abundance determination being biased by residual 
systematic errors. The following spectroscopic indicators were 
utilised in the quantitative analysis: 

• Teff: all available H and He lines, and multiple independent 
ionization equilibria; confirmation via spectral energy distribu- 
tions (SEDs); 

• log g: wings of all available hydrogen Unes and multiple ioni- 
zation equilibria; confirmation via Hipparcos distances; 

• ^: several elements with spectral lines of different strength 
enforcing no correlation between e(X) and the strength of the 
lines (equivalent to e(X) being independent of W,\)', 

• V sin / and ^: metal line profiles; 

• s{X): a comprehensive set of metal lines. 

The parameter determination started with the analysis of the 
hydrogen and helium lines. When a good simultaneous fit to 
most H and He lines was achieved - Tefr and \ogg were then 
typically constrained to better than ~5% and 0. 1-0.2 dex, respec- 
tively, for this high-quality set of stellar spectra - the procedure 
commenced to consider lines of other elements. Ionization equi- 
libria, i.e. the requirement that lines from different ions of an 
element have to indicate the same chemical abundance, facil- 
itated a fine-tuning of the previously derived parameters. The 
selection of ionization equilibria to be analysed depends primar- 
ily on Teff. Table|4]summarises the spectroscopic indicators em- 
ployed for the Teff and log g determination of each sample star, 
sorted by temperature. Elements that show lines of three ion- 
ization stages simultaneously in the spectrum are most valuable 
as they can in principle constrain both T^-s and \ogg at once. 
Examples are C ii/iii/iv or Si ii/iii/i\|_| in early B-type stars. When 
lines from only two ionization stages of an element are present in 
the spectra, then more indicators are required for the parameter 



' Note tliat we analyse Si iii/iv lines only because the model atom em- 
plo yed here underesti mates silicon abundances derived from Si ii lines, 



plo yeu nere unaeresti mates silicon abur 

Spectrum Plotting and Analysing Suite, Spas iHirschll2009h . see lSimon-DiazI l l201Cfl for a discussion 
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determination. Usually in the literature, a few hydrogen lines are 
analysed and one ionization equilibrium is established. On the 
other hand, we try to establish typically 4-5 independent ioniza- 
tion equilibria in addition to fitting all available hydrogen lines, 
which is unprecedented. Finally, the resulting model fluxes were 
compared with the observed SEDs. While such a high accuracy 
in the Teff-determination as with our spectroscopic approach can 
not be achieved by SED fitting alone, it provides a valuable con- 
sistency check, which can also facilitate the detection of other- 
wise unrecognised cooler companions. 

Microturbulence was for a long time an ad-hoc fit parame- 
ter that was employed to remove correlations of abundance with 
Wx, and often a different ^ was adopted for different elements. 
Only recently, a physical explanation for the phenomenon of mi- 
croturbulence in hot stars was suggested, likely being a conse- 
quence of subsurface convection (Cantiell o et al.ll2()()9l) . similar 
to the case of solar-type stars. Microturbulence needs to be con- 
strained simultaneously along with T^s and \ogg, such that its 
determination was a crucial part of our iteration procedure . The 
reason for this is that an inappropriately chosen microturbulence 
may lead to substantial shifts in Teff from the ionization equi- 
libria analysis, see Fig. 5 in iNieva & Przvbillal (1201 Obi) , which 
may remain unnoticed in the usual approach of using a minimum 
set of indicators for the parameter determination. Consequently, 
several elements were analysed for deriving the micro turbulent 
velocity. Our primary indicators were the Si, O and C lines, but 
the results were checked for consistency with the lines from the 
other elements as well. 

The detailed analysis of line profiles shows that rotational 
broadening alone is often not sufficient to explain the ob- 
served line shapes in hot stars. Agreement can be achieved 
when introducing a radial-tangential anisotropic macroturbu- 
lence (lGravl l2005'. p. 433f.) as additional broadening agent, see 
e.g. Fig. 11 of Przvbilla et al. (2006a). Consideration of macro- 
turbulence is therefore essential for meaningful line-fits using 
;^^-minimisation. A physical explanation of macroturbulence in 
hot stars was also only recently sugg ested, hkely being a collec- 
tive effect of (non-radial) pulsations ( Aerts et aDl2009 '). 

Usually in stellar analyses, once the stellar parameters are 
fixed one commences with the abundance determination, treat- 
ing this as an essentially independent step. In our approach the 
abundance and stellar parameter determination are tightly re- 
lated because of the use of ionization equilibria. In consequence, 
only few species (those not appearing as ionization equilibria 
in Table |4|i are left to finalise the analysis. Another difference 
to typical literature studies is the large number of spectral lines 
evaluated by us per species, and the consistency achieved from 
the different ionization stages of the various elements. All the 
various improvements in observations, modelling and analysis 
methodology facilitated analyses at much higher precision and 
accuracy to be achieved than possible in standard works. The 
quality of the analyses could be retained over a large param- 
eter space, spanning nearly 20 000 K in T^^ and ranging from 
close to the zero-age main sequence (ZAMS) to the giant stage. 
Consequently, an excellent match of the computed and the ob- 
served spectra is achieved globally and in the details, see the end 
of Sect.|5]for a discussion. 



5. Results 

5.1. Stellar parameters 

Table |5] summarises the stellar parameters derived from the 
quantitative spectral analysis. This includes the atmospheric pa- 



Table 5. Stellar parameters of the program stars" . 



# 


HD 


Teff 

lO'K 


(cgs) 


^ 


vsin i 
kms"' 


<r E(B - V) Vo 
mag mag 


Mo 


4pec 

pc 


^HIP 

pc 


1 


36591 


27.0 
± 0.3 


4.12 
0.05 


3 
1 


12 

1 




0.06 5.16 
0.00 0.01 


12.3 
0.3 


408 
26 




2 


61068 


26.3 
± 0.3 


4.15 
0.05 


3 

1 


14 

2 


20 

1 


0.07 5.49 
0.01 0.03 


11.6 
0.3 


434 
28 




3 


63922 


31.2 
+ 0.3 


3.95 
0.05 


8 

1 


29 

4 


37 
8 


0.09 3.82 
0.01 0.02 


18.9 
0.5 


389 

25 




4 


74575 


22.9 
± 0.3 


3.60 
0.05 


5 
1 


11 

2 


20 

1 


0.05 3.54 
0.00 0.01 


11.7 
1.2 


301 

24 


270 
10 


5 


122980 


20.8 
± 0.3 


4.22 
0.05 


3 

1 


18 

1 




0.01 4.32 
0.01 0.02 


7.5 
0.2 


150 
10 


156 
5 


6 


149438 


32.0 
± 0.3 


4.30 
0.05 


5 
1 


4 

1 


4 

1 


0.03 2.73 
0.01 0.02 


15.8 
0.7 


143 
9 


145 
11 


7 


886 


22.0 
± 04 


3.95 
0.05 


2 
1 


9 

2 


8 

2 


0.00 2.83 
0.01 0.04 


9.2 
0.3 


120 
8 


120 
8 


8 


29248 


22.0 
+0.25 


3.85 
0.05 


6 

1 


26 

2 


15 

5 


0.01 3.90 
0.01 0.04 


9.5 
0.3 


225 
15 


207 
8 


9 


35299 


23.5 
± 0.3 


4.20 
0.05 



1 


8 

1 




0.02 5.64 
0.01 0.02 


9.2 
0.3 


344 

22 


269 

24 


10 


35708 


20.7 
± 0.2 


4.15 
0.07 


2 
1 


25 
2 


17 

5 


0.06 4.69 
0.01 0.02 


7.6 
0.2 


192 
17 


192 
8 


11 


36512 


33.4 
± 0.2 


4.30 
0.05 


4 

1 


20 

2 


10 

5 


0.02 4.54 
0.01 0.03 


18.0 
0.7 


366 

24 




12 


36822 


30.0 
± 0.3 


4.05 
0.10 


8 
1 


28 

2 


18 

5 


0.11 4.06 
0.01 0.04 


16.2 
1.1 


338 
46 


333 
28 


13 


36960 


29.0 
+ 0.3 


4.10 
0.07 


4 

1 


28 

3 


20 

7 


0.02 4.74 
0.01 0.03 


144 
0.6 


392 
36 




14 205021 


27.0 
±0.45 


3.95 
0.05 


4 

1 


28 
3 


20 

7 


0.03 3.13 
0.01 0.03 


13.3 
0.5 


201 
13 


210 
13 


15 209008 


15.8 
± 0.2 


3.75 
0.05 


4 
1 


15 

3 


10 
3 


0.04 5.87 
0.01 0.04 


5.8 
0.2 


372 
25 




16 216916 


23.0 
± 0.2 


3.95 
0.05 




1 


12 

1 




0.08 5.33 
0.01 0.03 


9.8 
0.3 


405 
26 




17 


3360 


20.75 
± 0.2 


3.80 
0.05 


2 
1 


20 

2 


12 

5 


0.02 3.61 
0.01 0.03 


8.9 
0.3 


191 
12 


182 

5 


18 


16582 


21.25 
± 04 


3.80 
0.05 


2 
1 


15 

2 


10 

5 


0.00 4.07 
0.01 0.03 


9.3 
0.3 


241 
16 


199 
6 


19 


34816 


30.4 
± 0.3 


4.30 
0.05 


4 

1 


30 

2 


20 

7 


0.00 4.29 
0.02 0.05 


144 
04 


264 
18 


261 
16 


20 


160762 


17.5 
± 0.2 


3.80 
0.05 


1 
1 


6 

1 




0.00 3.80 
0.00 0.01 


6.7 
0.2 


157 
10 


139 

3 



Notes. *°* Tcff and log g are expected to vary over a pulsation cycle in 
the variable stars (cf. Table[T), see Sect. l5.1l for further discussion. 



rameters effective temperature Teff , surface gravity log g, micro- 
turbulence ^, projected rotational velocity v sin / and macrotur- 
bulent velocity (. Additional quantities include the computed 
colour excess E(B - V) of the sample stars, their de-reddened 
apparent magnitude Vo, evolutionary masses Mgy, spectroscopic 
distances lispec and Hipparcos distances duip. 

Effective temperatures are constrained to 1-2% and surface 
gravities to less than 15% (Icr uncertainties). It is unlikely that 
these values are subject to larger residual systematicsQ as they 
are constrained by the simultaneous match of many independent 



'" This is strictly valid only for the time of observation, see Sect. l2.2l 
The v ariable stars are expe cted to show (correlated) changes of T^g and 
log^, dde Ridder et alj200 2: Catanzaro & Leone 2008), which may ex- 
ceed the given uncertainties in Table |5] Average atmospheric param- 
eters as derived from the analysis of time-series observations may be 
more appropriate in this context, however our approach recovers the 
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Fig. 4. The sample stars in the Teff-log g plane. Icr-uncertainties 
are shown. Numbering according to Table [1] Overlaid are evo- 
lution tracks for non-rotating stars of metallici ty Z = 0.02 from 
[Meynet & Maeder (2003, full lines) and Sch aller et all (Il992i 
dotted lines). See Sect. lS.ll for a discussion. 



indicators. Even changes of the underlying physica l models, like 
a use of hydrostatic non-LTE model atmospheres (lNievall2007l) 
or hydrodynamic non-LTE model atmospheres (NSll), have 
been shown to have only small effects (i.e. agreement of stellar 
parameters and elemental abundances is obtained with the dif- 
ferent models, within the stated statistical Icr-uncertainties and 
without systematic trends). Overall, this is a major improvement 
over other literature studies, where the uncertainties can amount 
to ~5-10% for effective temperature and ~25% for surface grav- 
ity. In consequence, all quantities depending on temperature and 
gravity (e.g. reddening, evolutionary masses, spectroscopic dis- 
tances) also show reduced uncertainties. A slight degeneracy of 
line profile variations to simultaneous changes of v sin / and ( 
prevents one to achieve very tight constraints on these quantities 
individually. 

A noteworthy result is the finding that one value of micro- 
turbulent velocity is derived from the different chemical species, 
after several iterations in all variables of the spectral fitting pro- 
cedure. This was often not the case in previous studies, likely 
being a consequence of adopting ill-chosen atmospheric param- 
eters or of shortcomings in the modelling, e.g. the assumption 
of LTE or the use of limited sets of atomic data for non-LTE 
line-formation calculations. 

A comparison of the resulting model fluxes with the ob- 
served spectral energy distributions (where available) shows 
good agreement, see Fig. 4 (available online only). Our val- 
ues for effective temperature derived via multiple independent 
ionization equilibria are thus further verifiedF^ Th e photometric 
data w as converted into fluxes using zer opoints of Bessell et alj 
(Il998h for Johnson photometry and of iHeber et"alJ (2002) for 
the 2MASS photometry. All observed fluxes were de-reddened 
using an interstellar reddening law according to C ardelli et al.l 
(119891) . adopting colour excesses E{B - V) as indicated in Fig. 4 
and a ratio of total-to-selective extinction Ry -AyjEiB -V)- 
3.1. Note that the E(B - V) values used in Fig. 4 may differ 



time-invariant quantities like elemental abundances, which is the main 
topic of the present work. 

'* The opposite approach, using the SEDs as an independent Tcff- 
indicator is of limited value because of a much lower sensitivity of the 
method. 




400 



Fig. 5. Results of the distance determination. Upper panel: 
comparison of spectroscopic and Hipparcos distances. Icr- 
uncertainties are shown. Numbering according to Table [1] The 
1:1 relation is indicated by the dashed line. Lower panel: per- 
cent difference of the two distance indicators. Open circles mark 
the data if corrected for Lutz-KeUcer bias. See Sect. 15.11 for a 
discussion. 



slightly (on average by ~0.02 mag, which is within the mutual 
uncertainties) from those in Table|5] which were calculated from 
the difference between the observed and the computed Atlas9 
model colour As lUE spectrophotometry was not available for 
all the sample stars, we decided to use the homogeneously de- 
rived E(B - V) data in Table|5]to determine Vq- 

Evolutionary masses of the sample stars were determined by 
comparison of the objects' positions in an effective temperature 
Tes VS. s urface gravity log g diagra m with stellar evolution mod- 
els from iMevnet & Maedej (120031) . see Fig. |4] Tracks for non- 
rotating stars with 'solar' metallicity Z = 0.02 were adopted. As 
several stars are less massiv e than the lower ma ss limit of this 
grid, additional tracks from Schaller et al.l (Il992h were adopted, 
which however show a small offset with respect to the more 
modem data, see the 9 Mq models. For consistency, we applied 
corresponding offsets to the lower-mass objects in order to derive 
homogeneous Mev values in Table |5] Note that some systematic 
offsets will result due to difference between the model (Z - 0.02) 
and our derived metallicities (Z = 0.014). One consequence will 
be a shift of the zero-age main sequence towards higher gravi- 
ties, such that the high-gravity objects 6, 11 and 19 will fall on 
the ZAMS. We neglect this in the following, as the effects on the 
further discussion are small. 

Once these parameters are constrained, it is possible to de- 
termine spectroscopic distances dspec of the sample stars using a 
formula bv lRamspeck et al.l (1200 lb 



*spec 



7.11 X 10"* VA^ev//vlO("-4^»-l°S«' [pc] , 



(1) 



where Mev is expressed in units of Mq, the Eddington flux at 
the effective wavelength of the V filter Hy in ergcm"^ s"' Hz ' 
(derived here from the Atlas9 models), Vo in mag and log g in 
cgs u nits. The formula uses a flux calibration of Vega accord- 
ing to lHeber et al.l (11984) . The most crucial input quantity in the 
distance determination is the surface gravity. This opens up a 
possibility to independently verify our log g determination via 
comparison of the spectroscopic with Hipparcos distances as de- 
rived from parallaxes n from th e new reduction of the Hipparcos 
catalogue (Ivan Leeuwenll2007h . see Fig. \5\ Lutz-Kelker correc- 
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Table 6. Metal abundances for the sample stars. 



# 


HD 


C 


N 





Ne 


Mg 


Si 


Fe 


1 


36591 


8.33±0.08 (30) 


7.75+0.09(61) 


8.75+0.11(53) 


8.09+0.08(21) 


7.58+0.10(6) 


7.50+0.04(5) 


7.48+0.11(32) 


2 


61068 


8.27±0.07 (23) 


8.00+0.12(61)" 


8.76+0.09 (49) 


8.07+0.11(17) 


7.56+0.03 (3) 


7.53+0.06(5) 


7.51+0.11(28) 


3 


63922 


8.33±0.07(19) 


7.77+0.08 (23) 


8.79+0.10(39) 


8.07+0.07 ( 8) 


7.60+0.01 (2) 


7.49+0.12* 


7.51+0.08 ( 9) 


4 


74575 


8.37±0.10(19) 


7.92+0.10(56)" 


8.79+0.08 (45) 


8.05+0.08(12) 


7.51+0.10(6) 


7.52+0.12* 


7.51+0.09(27) 


5 


122980 


8.32±0.09 (22) 


7.76+0.08 (47) 


8.72+0.06 (52) 


8.07+0.07(14) 


7.50+0.05 (4) 


7.25+0.04(4)" 


7.44+0.11(27) 


6 


149438 


8.30±0.12(32) 


8.16+0.12(73)" 


8.77+0.08 (49) 


8.14+0.07(18) 


7.62+0.03 (3) 


7.52+0.06(2) 


7.54+0.09(21) 


7 


886 


8.37±0.08(17) 


7.76+0.07 (40) 


8.73+0.11(52) 


8.11+0.08(13) 


7.61+0.05 (4) 


7.38+0.03 (5) 


7.51+0.07(30) 


8 


29248 


8.29±0.13(15) 


7.93+0.09(41)" 


8.78+0.09 (47) 


8.07+0.07(12) 


7.55+0.08(3) 


7.54+0.06 (5) 


7.52+0.08 (25) 


9 


35299 


8.35±0.09(16) 


7.82+0.08 (40) 


8.84+0.09 (52) 


8.07+0.10(14) 


7.53+0.06 (4) 


7.56+0.05 (5) 


7.53+0.10(28) 


10 


35708 


8.30±0.09(15) 


8.22+0.07(38)" 


8.82+0.11(45) 


8.06+0.09(12) 


7.65+0.02 (4) 


7.51+0.03(5) 


7.58+0.06 (24) 


11 


36512 


8.35+0.14(19) 


7.79+0.11(22) 


8.75+0.09 (39) 


8.11+0.07(11) 


7.50(1) 


7.54+0.07 (2) 


7.53+0.03 ( 3) 


12 


36822 


8.28±0.14(22) 


7.92+0.10(31)" 


8.68+0.10(39) 


8.06+0.09(14) 


7.54(1) 


7.56+0.07 (2) 


7.52+0.04 ( 9) 


13 


36960 


8.35±0.09 (20) 


7.72+0.11(36) 


8.67+0.08(41) 


8.13+0.11(13) 


7.62(1) 


7.56+0.07 (2) 


7.48+0.09(13) 


14 


205021 


8.24±0.06(10) 


8.11+0.11(33)" 


8.64+0.13(44) 


8.17+0.11(10) 


7.53(1) 


7.50+0.09 (2) 


7.55+0.10(20) 


15 


209008 


8.33±0.09(10) 


7.80+0.11(18) 


8.80+0.11(21) 


8.02+0.11(14) 


7.51+0.07(4) 


7.42+0.04(4) 


7.53+0.08 (26) 


16 


216916 


8.32±0.07(17) 


7.78+0.10(40) 


8.78+0.08 (47) 


8.10+0.11(14) 


7.54+0.06(5) 


7.51+0.05(5) 


7.50+0.08 (21) 


17 


3360 


8.31±0.08(14) 


8.23+0.07(37)" 


8.80+0.08 (38) 


8.11+0.08(12) 


7.56+0.04(3) 


7.60+0.07 (5) 


7.55+0.07(19) 


18 


16582 


8.21±0.09(16) 


8.23+0.08(39)" 


8.79+0.07 (44) 


8.05+0.09(12) 


7.54+0.05 (4) 


7.50+0.05 (5) 


7.56+0.10(27) 


19 


34816 


8.38±0.05(10) 


7.81+0.15(29) 


8.71+0.09(35) 


8.18+0.05 ( 8) 


7.60(1) 


7.54+0.06 (5) 


7.54+0.07 ( 8) 


20 


160762 


8.40±0.07(13) 


7.89+0.12(39) 


8.80+0.09 (29) 


8.05+0.07(13) 


7.56+0.06 (4) 


7.51+0.05(4) 


7.51+0.08(22) 



Notes. Uncertainties represent the line-to-line scatter The number of lines analysed per element/star is given in brackets. *"' Excluded from the 
discussion in Table[7]and Fig.[T2] See Sect.Ofor details. <*> Adopted from PNB08. 



tions dLutz&Kelkeiil 19731; ISmithll2003"h have not been adopted 
to refine the Hipparcos parallaxes, as they should not be applied 
to measurements of individual stars (Ivan Leeuwenll2007l p. 87). 
Their potential impact is nevertheless rather small for the present 
sample stars, as visualised in Fig. |5] Overall, good agreement 
of the spectroscopic and Hipparcos distances is found within 
the unceitainties, except for the objects 9, 18 and 20, which, 
however, are still within the 3cr limits. An apparent systematic 
trend of increasing difference d^pec - dmp with increasing dis- 
tance (lower panel of Fig. |5j becomes marginal when only one 
object, #9, is disregarded, i.e. the regression line is then compat- 
ible with slope and offset zero. 



5.2. Chemical abundances 

Metal abundances of the sample stars are summarised in Table |6] 
where also the standard deviation from the line-to-line scatter 
and the number of analysed lines are indicated. The values are 
averages over all lines of a given species, giving each line in the 
different ions equal weight. A precision of the results of better 
than 25% is indicated typically. The individual line abundances 
are listed in Table 7 (available online only), where further details 
on the line formation calculations are also given: central wave- 
lengths A of the spectral lines, excitation energy of the lower 
level x, oscillator strengths log gf and the accuracy and source 
of the oscillator strengths. 

We deviated from the above procedure for the derivation of 
the helium abundances. Helium is the second most abundant ele- 
ment and therefore cannot be treated as a trace element: changes 
in its abundance modify the mean molecular weight of the at- 
mospheric plasma and thus can affect the atmospheric structure 
noticeably. However, the comparison of the observed with the 
computed spectra showed that the sample stars are indeed well- 
described by a (protosolar) e(He) = 10.99+0.05 (except for ob- 
ject 4, with a - still compatible - e(He) - 10.94+0.05). 



The metal abundances show a small scatter around the av- 
erage sample abundances, which is visualised later, in Fig. [12] 
The only peculiarities are enhanced nitrogen abundances in sev- 
eral sample stars, which can be understood in the framework of 
mix ing of CN-cyc led material into the atmospheric layers (see 
e.g. lPrzvbilla et al. ' 2010). The silicon abundance in object 5 is 
also conspicuously low, which may be an indicator for the on- 
set of chemical differentiation in the otherwise normal star. We 
decided to discard these peculiarities from further analysis. 

The most relevant sources of systematic uncertainties in the 
chemical abundances that arise from the spectral modelling, 
analysis and observations in our approach have akeady been dis- 
cussed in a series of papers. We refer to the work summarised 
in Table [3] for estimates of effects and systematics due to un- 
certainties in the atomic input data. Systematics due to basic 
model atmosphere input physics were investigated in particular 
by NP07 and NSl 1, while systematics due to d etails of analysis 
strategi es were discussed bv iNieva(2007) and Nieva & PrzvbiUal 
dloTOallbl) . 

Since our observational material leaves little room for obser- 
vational bias we concentrate here only in quantifying the sys- 
tematic effects on chemical abundances as introduced by un- 
certainties in the stellar parameters. We exemplify the effects 
of independent effective temperature, surface gravity and mi- 
croturbulence variations on oxygen and silicon abundances (the 
most sensitive among the elements in this parameter range) for 
the star #7 (HD886, yPeg) in Table 8. Two cases are investi- 
gated, for our and also for typical values in the literature. Note 
that this case represents our largest relative uncertainty in Tetj, 
one of the largest in A^/^, and typical uncertainty in log g, such 
that the example constrains the maximum systematic effects ex- 
pected for our sample stars. As correlations exist between pa- 
rameters - e.g. a higher Teff implies a higher \ogg in analyses 
- the true systematic uncertainties are hard to quantify in de- 
tail. Based on the given example and previous experience we 
estimate systematic uncertainties of the elemental abundances in 
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our sample stars, accounting for all factors, to be about 0.15 dex 
at maximum, with the majority being accurate to within better 
than about 0.10 dex. In comparison the systematics due to stellar 
parameter variations for typical uncertainties from the literature 
are much higher (Table 8), amounting to a factor ~2 in individ- 
ual cases. Usually, uncertainties in Teff are the most critical, but 
note the high sensitivity of the (rather strong) silicon lines to 
micro turbulence variations. 

A test for residual systematics can be made by searching 
for trends among the elemental abundances as a function of at- 
mospheric parameters. Figure 7 (available online only) displays 
metal abundances from the present work and from nine addi- 
tional Orion OB 1 stars of NS 1 1 (analysed in the same manner) 
as a function of Teff and logg. All data points cluster tightly 
around the average sample values and the respective Icr-error 
margins (see Table|7]i, except for the few cases discussed earlier, 
mostly for nitrogen. No significant trends either with T^ff or log g 
are found. 

Finally, we want to briefly comment on several potential 
sources of systematic uncertainties due to non-standard input 
physics in the context of model atmospheres. Potential weak 
stellar winds present have a negligible effect on the photo- 
spheric line spectrum, as the effects of the velocity field on the 
atmospheric stratification become noticeable only outside the 
line-formation region. Wind variability by larger amounts than 
the mass loss-rates of M<IQ^^Mq typical for B-type main- 
sequence stars has no effects on th e metal line spectrum in su- 
pergiants (e.g. lSchiller & Przvbilla'2008). Systematic effects on 
the abundance analysis due to the presence of magnetic fields 
are also not expected. Only three sarnple st ars have a confirmed 
magn etic field, ySCep (Henrichs et al."2000), f Cas (iNeiner et al.1 
l20()3h and rSco (Donati et al. 2006). To date, none of fliese 
shows observational evidence for abundance spots or vertical 
chemical stratification. 



5.3. Global spectrum synthesis 

Stellar parameters in Table |5] in combination with averaged 
chemical abundances in Table |6] were used to compute global 
synthetic spectra to visualize the quality of the analysis proce- 
dure. Overall, excellent agreement is found for all stars over 
the entire observed spectral regions. Examples for four stars 
with spectral types B3III (#15 HD 209008, 18 Peg), B2IV (#7 
HD 886, r Peg), Bl V (#1, HD 36591, HR 1861) and B0.2 V (#6 
HD 149438, t Sco) are shown in Figs. 8a-l le (available online). 

It is worth to notice that it is relatively simple to achieve 
reasonable fits of models to lower resolution spectra and data 
at lower S/N over limited wavelength regions, i.e. whenever the 
observational details tend to be washed out. High-resolution and 
high-S/N observations spanning a large wavelength range are 
much more challenging to be reproduced with synthetic spec- 
tra based on physical models at once. This is feasible only if the 
models match the global physics (i.e. the atmospheric structure) 
and details of individual features (i.e. the lines) equally well. 
We consider our success as a strong support for the absence of 
any significant systematics from our analysis. We thus comple- 
ment the probably most comprehensive benchmark test for stel- 
lar atm osphere modelling of OB stars to date bv Marcolino et all 
(l2009l at about the upper red-boundary of the present work). 

The locations of numerous spectral lines that are considered 
in our line-formation computations are indicated in the upper 
parts of the individual panels in Figs. 8a-l le in order to facilitate 
an evaluation of their presence/absence for a given set of stellar 



parametero- These include many more features (also from ad- 
ditional chemical species) than those analysed quantitatively in 
the present work (the latter are marked in the lower parts of the 
panels). We thereby show that even complex blend features like 
e.g. around Hei A412QA or around O ii/C iii /J4650A can be re- 
produced well. 

There are still some residual minor shortcomings. As can be 
seen from the comparison of observation with theory the con- 
tinuum normalisation could be improved globally. However, we 
have corrected for this by adjusting the continuum locally when 
analysing individual wavelength regions, such that this would be 
rather a cosmetic improvement. Our compromise to fit the five 
displayed Balmer lines. He to Ha, results in a slight mismatch 
in some of the Balmer line wings because of the imperfect nor- 
malisation, but the effect is well within the uncertainties of the 
log ^-determination. The cores of the Balmer lines for the hot 
star tSco (and just visible in HR 1861) are not matched well 
because non-LTE effects on the atmospheric structure are likely 
to affect the outer photosphere, i.e. the core-formation region 
for the strongest lines, gradually with increasing temperature. 
Moreover, Ha is expected to be influenced by the weak stellar 
wind, which is unaccounted for by our hydrostatic modelling. 
The unavailability of realistic broadening data is an issue for 
some lines, e.g. for Hei/13926 A. 

Several (high-excitation) spectral lines of the elements from 
Table [3] and some elements/ions are still not incorporated in 
our non-LTE spectral synthesis. On the other hand, we have 
included Al iii and S ii/iii in the calculations, adopting model 
atoms of Dufton et al. ( 1986) and Vranck en et al. (1996), respec- 
tively, and assuming solar abundances of lAsplund et al.l ( l2009l) . 
However, we do not consider these elements for the analysis as 
several shortcomings have been identified in the model atoms. 
In consequence, not all of the computed lines of these two el- 
ements give a good match to the observed spectra. A thorough 
verification and improvement of the model atoms for aluminium 
and sulphur is beyond the scope of the present paper, as is the 
inclusion of the missing lines from other chemical species. We 
will report on our efforts for achieving completion of the non- 
LTE spectrum synthesis in forthcoming papers. 



6. Present-day cosmic abundances 

The high degree of chemical homogeneity of the sample stars 
(Table |6]l encourages us to identify the average abundances 
with the long-sought cosmic abundance standard (CAS) for the 
present-day chemical composition of the cosmic matter in the 
solar neighbourhood. As statistical significance does matter for 
such a claim, we seek to compensate for the seven stars that 
were removed from our initial sample in order to avoid obser- 
vational biases. We therefore adopt abundance data for 9 addi- 
tional early B-type stars from the Ori OB 1 associatiorQ (NS 11), 
which were derived using the same models and analysis tech- 
niques as applied here. The stars meet the same observational se- 
lection criteria (i-iv in Sect. 12.1! ) as our core sample. The spectra 
were obtained with Fies on the 2.5m Nordic Optical Telescope 



'- More comprehensive line i dentifica ti ons f or early B-type 
stars are provided e.g. by Kili an et aU ( Il99lh for the optical 
blue spectral region, or by Gummersbach & Kaufer, online via 
http: //www. Isw.uni-heidelberg . de/cgi-bin/websynspec . cgi] 

'^^ Four other stars are in common with the present work: HD 35299, 
HD 36512, HD 36591 and HD 36960. The stellar parameters and abun- 
dances of NS 11 as derived from Fes spectra agree very well with the 
present data, providing another independent consistency check. 
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solar neighbourhood. Red histograms: present work, establishing the cosmic abundance standard. Black histograms: literature data. 
Photospheric and protosolar abundances from AGSS09 are also indicated, the bars representing the range spanned by the +lcr- 
uncertainties. See Sect. |6]for details. 



(La Palma), covering the wavele ngth range of 3700-7300 A at 
R^46 000 and S/N> 250 dSimon-Diaz 2010). 

Distribution functions for the individual elemental abun- 
dances in the star sample are displayed in Fig. [12] (red his- 
tograms). The data are normalised by the number of sample 
members, with the maximum value set to 1. The bin width is 
chosen as the standard deviation of the individual distributions. 
Note that only 20 stars are considered in the histograms for N 
- the atmospheres of nine stars are mixed with CN-processed 
material - and 28 for Si - one star is Si peculiar - as we are 
interested in the pristine abundances for constraining the CAS 
(see Table |6] for an identification of the data removed from the 
discussion here). 

Very tight distributions are found, described by a standard 
deviation of typically ~ 1 0% . Mean abundances together with the 
standard deviations (of the sample) are given in Table |7] which 
summarises also data on present-day abundances in the solar 
neighbourhood from other object classes, and the Sun. Resulting 
mass fractions for hydrogen (X), helium (Y) and the metals (Z) 
are indicated in Table [8] In addition to the metals investigated 
here - which cover the most abundant ones in the cosmos -, data 
for all other metals up to zinc was consi dered for cons training 
Z, using solar meteoritic values of Asplu nd et al.l (l2009l) . except 
for chlorine and argon, where abundances from the Orion nebula 
were adopted (Esteban et al. 2004). Any deviations of this aux- 
iliary data from the 'true' cosmic values will be absorbed by the 
error margins of Z due to their small contribution. 

This finding of homogeneity is in analogy to PNB08, but 
on a statistically much more significant basis and corrected for 
slight systematics in the iron abundances (see Sect. 13. 11 1. The 
same degree of chemical homogeneity is recovered as for the 
gas-phase of the diffuse ISM out to distances of 1.5 kpc from 
the Sun (Sofia & Mever 2001), see also Table|7] though different 
absolute abundance values for many elements are found because 
of depletion onto dust grains (see Sect. l7.3l l. 



Table 8. Mass fractions for hydrogen, helium and metals. 







Cosmic Standard 


Sun-^ 


- photospheric values 




B stars - this work 


GS98 


AGSS09 


CLSFBIO 


X 
Y 
Z 


0.710 

0.276 

0.014+0.002 


0.735 
0.248 
0.017 


0.7381 
0.2485 
0.0134 


0.7321 
0.2526 
0.0153 



However, the finding is at odds with practically all previ- 
ous work on early B-type stars in the solar neighbourhood. 
We concentrate on literature data from homogeneously anal- 
ysed samples of more than 10 stars for the comparison with 
the present work, as a comprehensive review of all available 
data is beyond the scope of this pape r. Therefore, abundances 
from the studies of lKilianI i 1 992l 1 1 994 . iGies & L ambertl (< 1 992|) 
- exclu ding bright g iants a nd sup er giants - -, |C unha & Lamb^ 
(1994), Daflon et al. (1999", 200 la'b", 2003"), Cunha et al. (200C 
Morel & Butler (2008) and Lyubimkov et al. ( 2004, 2005) are 
adopted, essentially applying the same distance cut as in our 
sample selection. These abundances were derived from high- 
resolution spectra using comparable non-LTE techniques as 
utilised here: either hybrid non-LTE or full non-LTE modelling 
under consideration of metal line blanketing. The only excep- 
tion are iron abundances, which were determined in LTE in these 
studies. The combined distribution functionofor the individual 
elemental abundances from the literature are also displayed in 
Fig. [T2l (black histograms). Much broader distributions are indi- 
cated, with typical standard deviations of about 0.2 dex. This is 
a factor ~5 larger than in the present work, despite our sample is 
a representative sub-set of the previously investigated stars, see 
Sect.lO 



'"* Note that many stars were subject to two or more of these indepen- 
dent investigations. No attempt is made to single out these cases: each 
analysis is considered with equal weight in the histograms. 
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Table 7. Chemical composition of different object classes in the solar neighbourhood. 



Cosmic Standard 



Orion nebula 



Young 



ISM 



Sun* 



Elem. 


B stars - this work" 


Gas 


Dust" 


F&G stars'' 


Gas 


Dust^ 


GS98 


AGSS09 


CLSFBIO 


He 


10.99±0.01 




10.988+0.003* 












10.93+0.01 




C 


8.33+0.04 


214+20 


8.37+0.03" 


~0 


8.55+0.10 


7.96+0.03^ 


123+23 


8.52+0.06 


8.43+0.05 


8.50+0.06 


N 


7.79+0.04 


62+6 


7.73+0.09* 






7.79+0.03^ 


0+7 


7.92+0.06 


7.83+0.05 


7.86+0.12 





8.76+0.05 


575+66 


8.65+0.03'- 


128+73 


8.65+0.15 


8.59+0.01'' 


186+67 


8.83+0.06 


8.69+0.05 


8.76+0.07 


Ne 


8.09+0.05 


123+14 


8.05+0.03'- 










8.08+0.06 


7.93+0.10 




Mg 


7.56+0.05 


36.3+4.2 


6.50:'- 


33.1+4.2: 


7.63+0.17 


6.17+0.02' 


34.8+4.2 


7.58+0.05 


7.60+0.04 




Si 


7.50+0.05 


31.6+3.6 


6.50+0.25'- 


28.4+4.3 


7.60+0.14 


6.35+0.05' 


29.4+3.6 


7.55+0.05 


7.51+0.03 




Fe 


7.52+0.03 


33.1+2.3 


6.0+0.3" 


32.1+2.5 


7.45+0.12 


5.41+0.04' 


32.9+2.3 


7.50+0.05 


7.50+0.04 


7.52+0.06 



Notes. *"' Including nine stars from Orion (NSII), in units of log(El/II) + 12/atoms per 10* H nuclei - computed fro m average star abun - 
dances (mean values over all individual lines per element, equal weight per line), the uncertainty is the standard deviation; '*' lEsteban et alj ( 120041) : 
*"'[Sim6n-Diaz & Stasinska ( 201 1); ''" difference between the cosmic standard and O rion nebula gas -phase abundances, in units of atoms per 10* H 
nuclei; *'' Sofia & Mevej ( 120011) ; ^^ value determined from strong-line transitions (S ofia et al.ll2011h . which is compatible with data from the anal- 
ysis of the [C ii] 158 ;um emission ( IDwek et al. 1997). Weak-line studies of C ii] /f2325 A indicate a higher gas-phase abundan ce £(C) =8.1 1+0.0 7 
dSofi a 2004), whic h corr esp onds to 84+28 ppm of carbon locked up in dust; '*' Mever et al. ( 1997), corrected accordingly to lJensen et al.l (l2007h ; 
*''* ICartledge et aU ( 120041) ; ^'^ ICartledge et al.l ( l2006h . The uncertainty in the ISM gas-phase abundances is the standard error of the mean; ^^' differ- 
ence b etween the cosmic standard an d ISM gas-ph ase abundances, in units of atoms per 10* H nuclei; **' photospheric values of lGrevesse & Sauvall 
( fl998l GS98), lAsDiund et alj ( l2009l AGSS09) and lCaffau et aLl ( l20ra CLSFGIO). 



If is exfremely difficult to trace the discrepancies in sfellar 
parameters and elemental abundances of individual stars from 
the various investigations, as differences exist at all levels: the 
qualify of the observational material, the methodologies of stel- 
lar parameter determination, the choice of analysed lines, the 
input atomic data, the computer codes and assumptions used for 
the modelling, among many other details that differ from study to 
study. A combination of several factors is most likely responsible 
for the discrepancies. We do not aim at resolving these discrep- 
ancies in detail case by case as little can be learned in terms of 
the objective of the present work. Moreover, the overall picture 
can actually be understood rather well from some basic consid- 
erations. 

Any abundance determination using a method with finite 
precision will yield an abundance distribution with a larger dis- 
persion than the true one. Broad distributions like those derived 
from the literature data can result from underlying tight distri- 
butions. Actually, the finding of such broad distributions is ex- 
pected, given that the statistical Icr-uncertainties in the literature 
data alone can amount fo 0.2 to 0.3 dex and systematic uncer- 
tainties result in shifts of the derived abundance distributions of 
the individual studies relative to each other (see Fig. 2 of PNB08 
for a visualisation). A full spread of the literature data over 1 dex 
and a shift in mean abundances is therefore hardly surprising. 

Finally, Fig. [12] allows also a comparison of the abundance 
distributions for the B-stars with the solar s tandar d to be made. 
Photospheric abundances of lAsplund ef al.l (1200 9*) are chosen as 
a representative for the type of data that is typically adopted in all 
kinds of astrophysical literature, and protosolar values from the 
same source as a representative for the bulk composition of the 
Sun (correcting for a ~0.04 dex depletion of the photosphere due 
to diffusion). Interestingly, similarities are found for some ele- 
ments and differences for others, which will be discussed further 
in Sect. 17.51 Note that the Icr-uncertainties of the solar abun- 
dances are about the same as the B-star sample standard devi- 
ations. Technically, the fact that a larger number of B-stars is 
considered for the determination of the CAS - in contrast to one 
star defining the solar standard -, means that the uncertainties 
of the CAS-values can be expressed via the respective standard 
error of the mean, which amounts to 0.01 dex for all metals stud- 



ied here. However, we prefer to assign the standard deviation of 
the sample as conservative error margins for the CAS. 

In our quest to reduce systematic eiTors to a minimum we 
cannot ignore possible bias introduced by other factors than stel- 
lar atmospheres alone. When considering an extended region 
like the solar neighbourhood (in our definition), effects from 
Galactic chemical evolution may also come into play. The pres- 
ence of Galactic abundance gradients implies a decrease of metal 
abundances with increasing Galactocentiic radius. To verify this, 
we checked for correlations of the stellar abundances with the 
spatial positions of the stars. An example for oxygen is shown in 
Fig. [13] No correlations are found, neither with Galactocentric 
distance nor with distance from the Galactic plane. We conclude 
that signatures of Galactic chemical evolution are insignificant 
on scales of ~500pc, providing no bias to the CAS on the level 
of precision achieved with our analysis methodology. 

In summary, it emerges from our previous discussion that the 
drastic reduction of many systematic uncertainties in our analy- 
sis procedure is the key for understanding the derived small dis- 
persion in the elemental abundances of the sample stars. For the 
first time the tiTie abundance distributions of the young stellar 
population in our Galactic vicinity are approximated, which ap- 
pear intrinsically tight. At the same time, the overall match of 
a large number of independent observational constraints and the 
successful passing of numerous tests for remaining biases puts 
confidence in the accuracy of the results. This allows an accurate 
and precise cosmic abundance standard to be established. We ex- 
pect that the true abundance distributions for our B-star sample 
will be in fact tighter than derived here because they are subject 
to analysis with a methodology of finite precision. Therefore, 
our study gives only a upper limit on the true degree of chem- 
ical homogeneity of the present-day cosmic matter in the solar 
neighbourhood. Given the current state of input physics for the 
models, it will be highly costly to improve the analysis inventory 
to a degree where much tighter constraints can be achieved. 

7. Implications 

In the following we want to investigate what implications the use 
of the cosmic abundance standard instead of the solar standard 
has on various astrophysical fields. We concentrate on the impact 
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Fig. 13. Spatial distribution of oxygen abundances in the sample 
stars, in analogy to Fig. [T] Open circles represent objects from 
NS 1 1 . The symbol size encodes the abundance according to the 
figure legend. 



of our sample data and the resulting CAS values for the evolu- 
tion of massive stars, for ISM science and for Galactic chemical 
evolution. Finally, we briefly comment on the origin of the Sun 
and its relation to its present Galactic neighbourhood. 

7. 1 . Stellar evolution: the Initial chemical composition 

The initial chemical composition has a profound influence on the 
structure and evolution of stars because of its effect on opacities 
and mean molecular weight. We have addressed the eff'ect of a 
metaUicity reduced from the so far canonical value Zq = 0.02 to 
ZcAS = 0.014 in the discussion of Fig. g] The shift of the ZAMS 
towards higher gravities will also be accompanied by a shift of 
the evolution tracks towards higher effective temperatures. But, 
would there be a significant effect if the CAS or modem solar 
abundances at about the same Z (see Table |8]i were used? In 
terms of the position of the tracks in the Hertzsprung-Russell 
diagram - probably not; however, in terms of observable tracers 
for rotational mixing, certainly yes. 

Energy production in massive stars is governed by the 
CNO cycles throughout most of their lifetime and the nuclear- 
processed material may reach their surface layers through ro- 
tational mixing alr eady during their main sequence phase (e.g. 
iMaeder & Mevneti.2000; .Heger & Langer.2000) , opening up a 
very powerful diagnostic to test models of stellar evolution. The 
changes of the CNO surface abundances reflect the actions of the 
dominating CN-cycle initially, following a well-defined nuclear 
path. This can be analytically approximated by a straight line 
in the diagnostic N/C-N/O-diagram (see Fig. fl4 |. with a slope 
define d solely by the initial CNO abundances dPrzvbilla et alj 
I2OIOI) . regardless of the mass, initial velocity or other model de- 
tails. 

The present analysis of our sample stars (and the additional 
nine objects from NSl 1) facilitates the predicted trend to be re- 
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Fig. 14. Observational constraints on mixing of CNO-burning 
products in massive stars. Mass ratios N/C over N/O are dis- 
played. Black dots: present data; black circles: 9 additional ob- 
jects from NSll. Data from previous non-LTE analyses are 
given as grey symbols ( green onli ne) - plus signs: Kilian ( 199:^ 
triangles: Gie s & Lambert (ri992[); dia monds: . Cunha & Lambert 
(1994), 
(2008); 



Dafl onetal.1 (1 19991 12001 
'Hunt er et al.1 (l2009l 






squares: iMorel et al 
The predicted nuclear 



crosses: 

paths assumin g initial CAS abundan ces (red star) and solar abun- 
dances (0) of lAsplundet al.l (l2009l) are indicated by the long- 
dashed and dotted lines, respectively. Statistical uncertainties of 
0.2 dex in each element - which are typical for previous litera- 
ture data - result in an error bars more than twice as large than 
those from the present data. 



covered for a statistically significant sample for the first time, 
see Fig. [14] Most of our objects cluster around the pristine CAS 
values, i.e. they are unmixed, while about 1/3 of the stars show a 
mixing signature of varying magnitude, following the predicted 
nuclear path with d(N/C)/d(N/0) = 4.6 (for initial CAS abun- 
dances) tightly. Stellar evolution models based e.g. on the solar 
values by AGSS09 would predict a different nuclear path (with 
slope ~3.0) despite a rather similar bulk metallicity. 

The large scatter found by previous non-LTE analyses of 
early B-type stars in the solar neighbourhood (many objects are 
in common with our sample) and additionally in three Galactic 
clusters is also likely a consequence of the lower accuracy and 
precision achieved in these studies, as argued in Sect. |6] Most 
data points are in fact consistent with the predictions because of 
the larger error bars, but they are of limited use for testing stellar 
evolution models stringently. 

Further consequences of the use of different individual abun- 
dances will be modified yields. As these are key input for 
Galactic chemical evolution models, they have to be determined 
as realistically as possible. 

7.2. Chemical homogeneity: early B-type stars vs. ISM 

From studies of interstellar absorption lines of the cold gas it 
is known for a long time that the local ISM out to 1.5 kpc 
from the Sun is che mically homogeneous, to the 10%-level 
JSofia & Meveij|200ll) . This can be understood as a natural con- 
sequence of turbulent mixing acting on all scales, which is 
due to the large density variations of the gas, generated by 
a complex interaction of many factors like momentum injec- 
tion by stellar winds and supernova shocks, magnetic fields 
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Fig. 15. Comparison of the stellar O abundance distribution as 
derived in the present work (red histogram) with gas-phase abun- 
dances alo ng different sightlines of the diffuse ISM (thin blue 
histogram. ICartledge et al.ll2004l) . in analogy to Fig. [12] Mean 
abundances are indicated. 



and self-gravity. Theoretical investigations of metal abundance 
fluctuations in the ISM have until recently been based on 
order of magnitude arguments involving characteristic spa- 
tial scales and timescales for vario us tur bulent mixing pro- 
cesses , see e.g. iRov & KunthI ( Il995h . and IScalo & Elmegreea 
( |2004|) for a review. More recently, sophisticated 3D (magneto- 
)hydrodynamic simulations of the local ISM at high resolution 
(e.g. de Avillez & Breitschwerdt 2007) support the view of an 
efficient mixing of metals at wide ranges of scales, driven by 
turbulence. 

The young stars in the solar neighbourhood are expected to 
follow the chemical characteristics of the matter from which they 
were formed. Indeed, the present study shows that, independent 
of the location of the sample stars in the solar neighbourhood 
(see Fig.[T3ll - whether they reside in OB associations or in the 
field -, and also independent of their mass (~6 to 20 Mq) and 
hence their life-time (~5xlO^ to 5x10^ yr, respectively), all stars 
show practically the same chemical composition. The fluctua- 
tions around the mean are ~10%, and probably less. This inde- 
pendent verification of the results from investigations of the ISM 
gas is achieved for the first time. 

The huge advantage of studying early-type stars is that the 
entire metal content can be determined using quantitative spec- 
troscopy, with no material hidden in an observationally inacces- 
sible reservoir like the dust-phase in the ISM. This opens up the 
possibility to determine the chemical composition of the dust in 
an indirect way, see the next sub-section. 

Moreover, our results put constraints on the injection and 
mixing timescales of metals in the local Galactic ISM, and there- 
fore on the hydrodynamics of the ISM. It appears that fresh nu- 
cleosynthesis products from supemovae and AGB (super-)winds 
or infall of pristine material onto the Galactic disk are unlikely 
to lead to a noticeable (at the level of the present abundance de- 
termination precision) local enrichment or depletion in a high- 
metallicity environment like the solar neighbourhood over a 
timescale of several 10^ yr Or, in other words, the interaction 
of hydrodynamic mixing on the one hand and viscosity and 
molecular diffusion on the other is highly efficient, such that the 
medium is homogenised quickly. 



7.3. Dust composition of the local ISM 

An important open question in our understanding of the ISM is 
the chemical composition of the dust particles. The amount of 
metals incorporated into dust cannot be derived directly from 
observations. Only an indirect determination is feasible, from 
the comparison of the ISM gas-phase abundances and a suitable 
reference that is unaffected by depletion onto dust grains. There 
is an ongoing debate which kind of objects provide the reference 
suited best for the co mparison: young B-typ e stars, young F&G- 
type stars, or the Sun (ISofia & MeveHl2001b . 

We argue here that the CAS values as established from the 
analysis of early B-type stars provide the long-sought reference, 
unprecedented in precision and accuracy. Most notable is that 
the same degree of chemical homogeneity is found for both, 
the CAS reference and the gas-phase abundances of the diffuse 
ISM (see Table|7]i. The abundance distributions are very similar, 
see Fig. [T5]for the example of oxygen (the gas-p hase abundance 
distribution is based on 37 diffuse sightlines of ICartledge et aU 
2004). From this it follows immediately that the dust-phase is 
chemically also rather homogeneous, what can be expected if 
mixing processes are highly efficient within the ISM (Sect. lTTSl i. 

We derive a similar chemical composition for the dust grains 
as PNB08, but at much higher statistical significance and with 
two exceptions (see Table|7|. The Fe abundance is higher due to 
the identification and removal of residual systematics in the line- 
formation computations for that elem ent (Sect. lTTT i. For carbon, 
a recent investigation by ISofia et aP (12011) raises doubts about 
the precision of weak-line analyses based on the C ii] /12325 A 
transition (see ISofial ( |2004|) for a discussion), possibly related 
to a systematically underestimated oscillator strength. We there- 
fore adopt a m ean abundance fro m five sightlines of the strong- 
line analysis of lSofia et al.l (1201 1[). which is compa tible with data 
from the [C ii]A15Sfim Une bv lDwek et alJ (Il997h . 

Overall, the results indicate a silicate/oxide-rich and rela- 
tively carbon-poor composition for the local ISM gas phase. In 
particular, the CAS provides sufficient oxygen to sustain the val- 
ues required by magnesium, silicon and iron to be locked up 
by vast majority in silicates (plus a small fraction in metal ox- 
ides) in the diffuse ISM. Using a reasonable dust composition 
(iDraind 12003) this amounts to about 140-150ppm of oxygen 
for the given abundances of the refractory elements in the dust, 
in good agreement with the observed value (Table |7]l, with some 
additional oxygen possibly incorporated in organic compound 
material. About 60% of the total carbon resides in the dust phase. 
Despite a higher abundance of carbon is found in the dust in ab- 
solute terms relative to PNB08, this falls still somewhat short 
of the demands of most dust models, see e.g. the discussions by 
Snow & Witt ( 1995) and Zubko et al. (2004). 

Finally, we want to comment on the cosmic abundance stan- 
dard in the context of gas- and dust-phase abundances in the 
Orion nebula. Further information for the general picture can 
be gained under the - not unlikely - assumption that the Orion 
giant molecular cloud formed out of material typical for the dif- 
fuse ISM, with subsequent chemical processing taking place in 
the cloud core. 

A compari son o f the CAS with gas-phase ab undances 
(lEsteban et al] l2004t ISimon-Diaz & Stasiriskal l201lh implies 
that the H ii region is devoid of carbonaceous dust. From ob- 
servations of the photodissociation region in the Orion nebula 
it is known that polycyclic aromatic hydrocarbons (PAHs) dis- 
appear as the gas becomes ionized (e.g. Tielens 2008). Our re- 
sults imply that photoevaporation affects all carbon-bearing dust 
particles in a similar manner, indicating that little carbon was 
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(2009, solar-type dwarfs and subgiants); diamonds: Spite et al. (2005, unmixed cool giants); plus signs: Israelian et al . (2004, un- 
evolved solar-type stars). Solar abundance ratios of Asplund et al. (2009) are also indicated (O). Right panels: black dots: unmixed 
objects from the present work, black circles: unmixed stars from NSll. Data from the literature like in the left panels. The CAS is 
also indicated (red star). Error bars (statistical 1 cr-uncertainties) typical for individual stars in the present study are shown. 



incorporated initially in graphite (the most stable form of car- 
bon under interste llar conditions), in line with the findings of 
lAmari et aP (Il990t) from laboratory studies of meteorites. Also, 
amorphous carbon dust grainsQ are either efficiently destroyed 
inside the ionized region, or they were a minority species ini- 
tially as well. On the other hand, there is only weak evidence for 
the destruction of silicate grains from the numbers in Table [Ti- 
the abundances of oxygen and of the refractory elements in the 
Orion nebula dust-phase are compatible with the ISM dust data 
within the (large) error bars. Unfortunately, a more direct verifi- 
cation of the carbon-poorness of dust within the H ii region via 
e.g. the extinction properties is complicated, as most of the ex- 
tinction tow ards Orion occurs in t he neutral medium surrounding 
the H II gas (Baldwin et al.'l99l'). 

The combined evidence from abundances in the ISM and 
in the Orion H ii region indicates that dust models considering 
silicates, PAHs, organic refractory material and possibly amor- 
phous carbon, but not graphite, should be investigated more 
closely. Models in analogy to th e COMP-NC-type or COMP- 
AC-type models of IZubko et al.l (12004.) look highly promising 
for future studies in view of abundance demands and the abil- 
ity to match other observational constraints like extinction and 
emission properties of the dust. We are confident that the tight 
observational constraints provided by the CAS will facilitate a 
better understanding of the nature of dust and grain structure to 
be developed. 

7.4. Galactochemical evolution: present-day abundances 

Nucleosynthesis in successive generations of stars has enriched 
the cosmic matter with heavy elements ever since the first 
Population III stars were bom. Studies of various objects like 
(Galactic and extragalactic) stars and H ii regions, or the ISM in 



damped Lya (DLAs) systems allow the cosmic enrichment his- 
tory to be traced and the specific production sites of individual 
elements to be constrained. The CAS provides valuable input 
for the comparison of models with observations, as it marks the 
present-day endpoint of galactochemical evolution, in particular 
for a typical spiral galaxy like the Milky Way. We put CAS val- 
ues into the context of the evolution of the five most important 
chemical species, the light elements CNO, magnesium as a typ- 
ical tracer of the a-process in core-collapse supernovae and iron 
as tracer of iron-peak nucleosynthesis in supernovae of type la. 

Observational constraints on the interlinked evolution of 
CNO are displayed in Fig. [161 Abundances from early B-type 
stars are compared to data from solar-type stars in the Galactic 
thin and thick disk, and in the halo. Typical statistical accura- 
cies for abundances in solar-type stars are ~0.05-0.10dex for 
LTE analyses, similar to the present work. Note, however, that 
systematic uncertainties due to non-LTE and 3D effects (and ad- 
ditionally due to the presence of magnetic field and stellar ac- 
tivity/chromospheres) can be substantial for solar-type stars, but 
are not understood comprehensively at presenO The compar- 
ison of literature data for the B-stars in the left panels and our 
data in the right panels shows once more which improvements 
in the precision and accuracy of the analyses were achieved. 

In terms of the investigation of the cosmic chemical evolu- 
tion the current focus of studies in the literature lies on the early 
phases at low metallicity. The interpretation of the data comes 
from the comparison with Galactic chemical evolution (GCE) 
models, which have to match the present-day composition as 
a boundary condition. It is therefore important for the entire 
modelling which reference values are used, solar or CAS abun- 
dances. In particular the differences in the C/O ratio are apprecia- 



" Amorphous carbon is considered the predominant grain material 
produced by C-type as ymptotic giant branch stars, the main source of 
carbonaceous dust (e.g. lWallerstein & KnaDdri998h . 



'' Approximate corrections for such effects may e.g. reproduce the 
upturn in the C/O ratio at low O as indicated by t he (non -LTE) data of 



iFabbian et al.l ( 12009. ) also in the data of Sp ite et aPdlOOST) . However, we 
prefer to display uncorrected values, as these are likely more realistic in 
view of more recent investigations (M. Spite, priv. comm.). 
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ble, also with respect to the majority of nearby solar-type stars, 
amounting to almost 50%. Taken at face value this difference in- 
dicates that the C/O enrichment of the ISM in the present-day 
solar neighbourhood occurred slower than at the Sun's place of 
birth (see Sect. 17.5b . No such conclusions could be drawn from 
the available literature data on early B-type stars so far 

Note that the Sun may be viewed as an extreme but still 
compatible case in terms of the distribution of the stars from 
the present sample in the logO-logC/O diagram, but the ab- 
solute values for the carbon abundance differ significantly. The 
solar and CAS data on N/0 are rather compatible on the other 
hand. A systematic investigation of nitrogen abundances in high- 
metalUcity solar-type stars would be desirable for further com- 
parisons. 

The occurrence of a-enhancement in the old stars of the halo 
and of the thick disk. Fig. [17] is well understood in terms of 
the different evolution timescales of supemovae of type la and 
II. The CAS Mg/Fe-ratio is compatible with the data from the 
Sun and from nearby solar-type stars, though it is somewhat 
low. Metal-poor objects among the solar-type stars are easily ex- 
plained by their large lifetimes. However, metal-rich stars are 
absent among the young stars from the present sample, while 
they are common among the solar- type stars. Finding a reason 
for this is not straightforward in terms of standard GCE (metal- 
licity is supposed to increase in time), but see the discussion in 
Sect. l7.5l related to radial migration. Again, data on early B-type 
stars from the literature were inconclusive in these terms. 

Overall, it is astonishing how different and at the same time 
how similar the young and old star populations in the solar 
neighbourhood are. It is for the first time that this is elaborated, 
as the lack of high precision and accuracy in many previous stud- 
ies of early B-type stars prevented any meaningful conclusions 
to be drawn. 

In addition, reference abundances are not only of interest in 
terms of the temporal evolution of elemental abundances, but 
also for the spatial distribution, in particular for the interpretation 
of Galactic abundance gradients. The current sample is not use- 
ful for deriving abundance gradients per se because of the small 
baseline spread in terms of Galactocentric radius (about 400 pc). 
However, the CAS provides a highly robust anchor point for the- 
oretical models at the solar circle, implying agreement with the 
solar standard for some elements, but also systematic shifts of 
various degrees for other elements. 

7.5. The Sun: an immigrant to its current neighbourhood 

In a strict sense, comparisons of the CAS and the solar chemical 
composition discussed in the previous sections are meaningful 



in terms of GCE only when there is a causal relationship of the 
4.56 Gyr old solar matter and the present-day material in its sur- 
roundings. In other words, the question is whether the matter that 
we see in the nearby new-born stars, as represented in particular 
by the early B-type stars, and in the local ISM has chemically 
evolved from the matter that was around at the birth of the Sun. 
The existence of such a relationship is a strong assumption for 
most previous and current GCE modelling^ 

Unfortunately, the birth-place of the Sun is unknown and 
tracing its orbit back is a highly complex task. Passages near 
other stars or molecular clouds, which are a stochastic process, 
and dynamical interactions with spiral arms prevent a straight- 
forward integration of the Sun's orbit in the Galactic potential 
backward in time from being successful. Improvements in our 
detailed understanding of the relevant processes affecting the 
motions of stars will certainly result from the Gaia mission in 
the future. But for the moment, we can rely only on a statistical 
approach, using the the oretical framework of Galactic dynam- 
ics and kinematics (e.g. lWielen et al.lll996l : [Sellwood & BinnevI 
2002) that predict that old stars like the Sun are able to migrate 
up to several kiloparsecs radially through the Galactic disk over 
their lifetime. It is therefore possible that the Sun was born away 
from the present solar circle. Chemical abundances can provide 
valuable additional constraints. 

The radial migration of stars within the Galaxy has been in- 
corporated i n Galactic chemical evolutio n calculations for the 
first time by ISchonrich & BinnevI (l2009l) . This kind of models 
constitutes a more reaUstic theoretical frame that could b e used 
to refine the previous estimations (e.g. Wielen et al. 1996) of the 
solar place of origin. It is important, though, to bear in mind 
that details in the model input like e.g. the adoption of different 
yields or abundance gradients along the Galactic disk may have 
a non-negligible effect in such estimations (R. Schonrich, priv. 
comm.). The required systematic studies are beyond the scope 
of the present paper, but we can provide some qualitative evalu- 
ation, which is illustrated by Table [9] 

In order to make a meaningful comparison of metal abun- 
dances in terms of Galactic chemical evolution, CAS values need 
to be compared to the bulk composition of the Sun (i.e. proto- 
solar value rT). corrected for the effects of GCE. The required 
abundance enrichments due to GCE are difficult to be quanti- 



" Technically, this is realised in GCE models by a division of the 
Galactic disk into concentric annuli that evolve independently from 
each other. The question in our context is, whether the Sun was bom 
at a similar Galactocentric radius as it is observed today. 

'^ Over the lifetime of the Sun, the combined effects of thermal dif- 
fusion, gravitational settling and radiative acceleration have lead to a 
build-up of helium and metal abundances below the convection zone. 
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Table 9. Chemical tagging of the Sun's place of birth. 



Element 



Protosun 



Protosun, GCE corrected" 



CAS 



AGSS09 CLSFBIO 



AGSS09 



CLSFBIO 



d£(El.)/di?g 
dexkpc"' 



CAS+d£(El.)/dj;g 
.Rg = 6kpc /?g = 5kpc 



C 

N 

O 

Mg 

Si 

Fe 



8.47+0.05 
7.87±0.05 
8.73±0.05 
7.64±0.04 
7.55±0.04 
7.54±0.04 



8.54±0.06 
7.90±0.12 
8.80±0.07 



7.56±0.06 



8.53+0.05 
7.95±0.05 
8.77±0.05 
7.68±0.04 
7.63±0.04 
7.68±0.04 



8.60±0.06 
8.01±0.12 
8.84±0.07 



7.70±0.06 



8.33+0.04 
7.79+0.04 
8.76+0.05 
7.56+0.05 
7.50+0.05 
7.52+0.03 



-0.103+0.018* 
-0.085+0.020^ 
-0.035'''' 
-0.039'' 
-0.045'' 
-0.052'' 



8.54+0.05 
7.96+0.05 
8.83+0.05 
7.64+0.05 
7.59+0.05 
7.62+0.03 



8.64+0.05 
8.05+0.05 
8.87+0.05 
7.68+0.05 
7.64+0.05 
7.68+0.03 



Notes , '"' applying values from Table 5 of AGSS09, based on G CE models of IChiappini et alj (l2003h : <*' lEsteban et all (l2005h : ^'^ ICarigi et al.l 
| |2005|) : ''^' ICescutti et al.l ( r2007h . based on Cepheid observations of'Andrievskv et all ( 120041 and references therein); *'"' a slightly steeper - though 
compatible - gradient, by -0.044+0.010 dex kpc"' , is given by £axigLet.al:. ( 2005h . 



fied, as they depend on many model details like e.g. the adopted 
star-formation history and yields. Crucial is that enrichment is 
expected to occur, not depletion. In view of the previous discus- 
sion in this section and in Sect. 17.41 it is therefore in fact the 
similarities of the solar values and the present-day CAS, which 
are astonishing, and less the differences. 

Let us elaborate the argument in more detail: assume that 
the protosolar nebula would not have collapsed to the Protosun 
(with elemental abundances according to columns 2 and 3 of 
Table |9]l- Instead, the gas would have been enriched in metal 
content over the past 4.56 Gyrs as predicted by GCE models for 
the present solar vicinity, leading to the formation of the Sun at 
the present day, with abundances according to columns 4 and 
5. The Sun would appear significantly more metal rich than its 
surroundings, represented by the CAS (column 6). 

A different birth place of the Sun than around the solar circle 
could resolve this apparent contradiction. Higher abundance val- 
ues have been characteristic for the inner disk of the Milky Way 
for a long time over Galactic history. This does not only apply to 
individual metal abundances, the argument is further sustained 
by the presence of a higher C/O ratio in the Sun than in the CAS 
(independent of the reference of solar abundances). The inner 
disk shows a higher C/O ratio, which is supported both by an ob- 
serva tionally derived negative radi al C/O gradient (lEsteban et alJ 
12005 ) and some GCE models (e.g. lCarigiet al.ll2005h . 

In order to further constrain the origin of the Sun, we 
have to correct the CAS data for the effects of Galactic abun- 
dance gradiento (column 7 in Table |9]l, which were derived 
from a carefully analysed sample of H ii regions (lEsteban et alJ 
l2005t ICarigietal.' '2005") and the modelling of the Cepheid 
data of [A ndrievskv et al. (2004, and references therein) by 
ICescutti et al.l (1200 7). From our own experience, we consider 
these w ork s as source s of reliable data on this topic, seel Przvbillal 
( |2008|) and iFimsteinI ([2010) for a discussion. The results from 
the gradient-corrected CAS values at Galactocentric radii R„ of 
6 and 5 kpc are shown in columns 8 and 9 of Table |9] 

The comparison of the GCE-coiTected protosolar abun- 
dances with the gradient-corrected CAS values suggests that the 
birthplace of the Sun was located at a Rg between 5 to 6 kpc (de- 
pending on the solar standard reference, which introduces an- 
other uncertainty to the previous considerations). We conclude 
that the Sun - and probably many other nearby older and metal- 
rich stars - are immigrants to the present solar neighbourhood, 
supporting views that stellar radial migration is essential for un- 



such that the photospheric abundances are not representative for the 
bulk composition of the Sun. 

" Data are for the present-day, but GCE models suggest that the flat- 
tening of Galactic abundance gradi ents over the lifetime of the Sun are 
insignificant for our considerations (iMarcon-Uchida et al.ll201(]l) . 



derstanding Galactic evolution dSchonrich & BiniievI 120091) . It 
would be highly interesting though difficult to investigate the 
details of the Sun's migration to its current location close to the 
Galactic co-rotation r adius, which is so favourable for the exis- 
tence of life on Earth (iLeitch & Vasishtiiil998i) . 



8. Summary 

With the present work we have established a new benchmark 
for analyses of early B-type stars, demonstrating that practically 
the entire observed spectra can be reproduced reliably. It was 
shown that by the combined use of sophisticated models and of a 
thorough analysis methodology on high-quality spectra absolute 
stellar parameters and elemental abundances for early-type stars 
can be determined with a precision rivaling that of differential 
studies of solar-type stars. The successful simultaneous match of 
many independent observational indicators like the Balmer line 
wings, multiple ionization equilibria, SEDs and the agreement 
of spectroscopic and Hipparcos distances indicates that high ac- 
curacy was achieved at the same time, facilitated by the minimi- 
sation of systematic errors wherever they could be identified. 

Overall, our analysis methodology provides both the accu- 
racy and precision to use early B-type stars as versatile tools 
for astrophysics, besides our immediate objective here. Some 
conclusions from earlier applications get even strengthened in 
retrospect by the present work; e.g. on chemical tagging in or- 
der to determine population membership o f hypervelocity stars 
(iPrzvbillaetal .11200831 iTilUch et al.ll2009l) . the investigation of 
subtle abundance signatures to constrain supernova nucleosyn- 
thesis from binary supernova runaways (iPrzvbilla et al.ll2008bt 
llrrgang et al.l 120101) . or on coupling quantitative spectroscopy 
with ast eroseismic analysis of pulsating OB stars observed with 
CoRoT ( Briquet et al.l20 1 f). And, the methodology offers a high 
potential for future applications, that may facilitate many more 
facets of astrophysics to be studied at high confidence. 

The principal application of our novel analysis methodology 
here was on a larger sample of nearby apparently slowly -rota ting 
early B-type stars that were cleaned from peculiar objects. This 
confirmed that the young stellar component of the solar neigh- 
bourhood is chemically homogeneous to better than 10%, in ac- 
cordance to studies of the local ISM abundances - the material 
out of which the stars formed. This in turn allowed us to estab- 
lish a present-day cosmic abundance standard (CAS), which has 
the advantage of redundancy when compared to the solar stan- 
dard, as it is based on an entire sample of stars instead of one 
object alone. So far, information for the most abundant chemical 
species have been provided, with log s > 7.50. 

First implications of the existence of a present-day cosmic 
abundance standard were outlined here. The CAS represents the 
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recommended initial chemical composition for stellar structure 
and evolution calculations, in particular for short-lived massive 
stars. Observationally, nitrogen enrichment in about 1/3 of the 
sample objects indicates mixing of the surface layers with CN- 
processed material from the stellar core, with the abundance ra- 
tios for the light elements carbon, nitrogen and oxygen following 
tightly the predicted nuclear path of the CNO cycles. 

The CAS provides the to date most authoritative reference 
for constraining the chemical composition of the ISM dust- 
phase. A silicate-rich and relatively carbon-poor nature of the 
local ISM dust is inferred, challenging many contemporary dust 
models. In combination with the finding that carbonaceous dust 
is practically absent inside the Orion Hii region this implies 
that the dust in the ISM is in its majority composed of silicates, 
PAHs, organic refractory material and possibly amorphous car- 
bon, with only little carbon incorporated in graphite. 

The CAS provides tight reference points for anchoring mod- 
els of the chemical evolution of the Milky Way, constraining 
the present-day endpoint of Galactic nucleosynthesis in the solar 
neighbourhood in the course of the cosmic matter cycle. We dis- 
cussed in particular the examples of the light elements CNO, and 
the Mg/Fe ratio as a tracer for the c/Fe ratio. While the CAS in- 
dicates a surprisingly good agreement with most inferences from 
the solar abundances, a striking difference is found for the C/O 
ratio, amounting to almost 50% between cosmic and solar val- 
ues. Intriguingly, early B-type stars with super-solar metallic- 
ity are absent in the solar neighbourhood, while several nearby 
solar-type stars show super-solar metallicities. This can be inter- 
preted as a consequence of radial migration of stars within the 
Galactic disk. 

Finally, the comparison of the CAS with the solar standard in 
view of Galactic chemical evolution (GCE) shows that the highly 
successful use of the Sun as proxy for cosmic abundances is 
somewhat coincidental. Radial migration outward from its birth 
place in the inner disk at ~5-6 kpc Galactocentic distance (where 
higher metallicity values were reached earlier in cosmic history) 
over its lifetime to its current neighbourhood has compensated 
for the expected metal enrichment in the course of GCE. A tell- 
taling signature is left only in the C/O ratio. The present work 
provides the so far most stringent evidence in terms of chemical 
signature that the Sun - like many other solar-type stars - is an 
immigrant to its current Galactic neighbourhood. 
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